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1 INTRODUCTION 



Over the past two decades, spintronics (jZutic et all 120041 ) has emerged as one of the most 
vigorously pursued areas of condensed matter physics, materials science, and nanotechnol- 
ogy. A glimpse at the Oxford English Dictionary reveals the following attempt to define the 
field succinctly: "spintronics is a branch of physics concerned with the storage and transfer of 
information by means of electron spins in addition to electron charge as in conventional elec- 
tronics." The rise of spint ronics was ignited by basic research on ferro magnet /normal- metal 
multilayers in late 1980s (jMaekawa and Shinjol . 120021 ). as recognized by the Nobel Prize in 
Physics for 2007 being awarded to A. Fert and P. Griinberg for the discovery of giant mag- 
netoresistance (GMR). The GMR phenomenon also exemplifies one of the fastest transfers 
of basic research results in condensed matter physics into applications where in less then ten 
years since its discovery it has revolutionized information storage technologies by enabling 
100 times increase in hard disk storage capacity. 

In recent years frontiers of spintronics have been reshaped through several intertwined 
lines of research: (i) ferromagnetic metal devices w here the main theme is manipulation of 
magnetization via electric currents and vice versa ([Ralph and Stilesl . 120081 ): (ii) ferromag- 
netic semiconductors which, unlike metal ferromagnets, offer additional possibilities to ma- 
nipulate their magnetic ordering (such as Curie temperature, coercive fields, and ma gnetic 
dopants), but are still below optimal o perating temperature (Jungwirth et al. . 20061 ): (Hi) 
paramagnetic semiconductor spintronics (lAwschalom and Flatta 120071 ) largely focused on all- 



electr ical manipulation of spins via spin-orbit (SO) coupling effects in solids ([Fabian et al 
20071 ): and (iv) s pins in semiconduct ors as building blocks of futuristic solid-state-based quan- 



tum computers ([Hanson et al. . 2007 ). Unlike early non-coherent spint ronics phenomena (such 
as G MR), the major themes of the "second-generation" spintronics (lAwschalom and Flatte . 
20071 ) are moving toward the spin coherent realm where spin component persists in the direc- 
tion transverse to external or effective internal magnetic fields. Recent experiments exploring 
such phenomena include: spin-transfer torque where spin current of large enough density 
injected into a ferromagnetic layer either switches its magnetization from one static con- 
figuration t o another or generates a dynamical situation with steady-state precessing mag- 
netization ([Ralph and Stilesl . 120081 1: spin pumping as the "inverse" effect of spin-transfer 
torque in which precessing magnetization of a ferromagnetic layer emits pure spin currents 



into adjacent normal metal layers in the ab sence of any bias voltage ([Saitoh et all 12006 



Cost ache et all 120061 ; iMoriyama et. all 120081 ) ; transport of coherent s pins (able to preces s 
in the external magnetic fields) across ~ 100 /xm thick silico n wafers ([Huang et all 120071 ): 
and the dire ct and inverse spin-Hall effects ( SHE) in bulk ([Kato et all 12004 qh_ and low- 



2006; 


Saitoh et al.. 


2006; 


Kimura et al.. 


2007) 



dimensional (IWunderlich et all 2005; S ih et all 120051 ) semiconductors and metals ([Valenzuela and Tinkhaml . 



SHE in both metals and semiconductors has 



been observed even at room temperature ([Stern et all 120061 ; iKimura et all 120071 )]. A closely 
relate d effort that per meates these subfields is the generation and detection of pure spin cur- 
rents ([Nagaosal . 120081 ) which do not transport any net charge, as illustrated in Fig. [TJ Their 
harnessing is expected to offer both new functionality and greatly reduced power dissipation^ 



The Joule heat losses induced by the current flow set the most important limits (jKeve sl l2005l) for conven- 
tional electronics, as well as for hybrid electronic-spintronic or purely spintronic devices envisioned to perform 
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Figure 1: The classification of spin I and charge I currents in metal and semiconductor 
spintronic systems corresponding to spatial propagation of spin-f and spin-j electronic wave 
packets carrying spin- resolved currents P and 1^: (a) conventional charge current I = I> + 
1^ 7^ is spin-unpolarized I s = — 1^) =0; (b) spin-polarized charge current / 7^ is 

accompanied also by spin current I s 7^ 0; and (c) pure spin current I s = ^(/ T - I 1 ) ^ 
arising when spin-| electrons move in one direction, while an equal number of spin-j electrons 
move in the opposite direction, so that total charge current is I = 0. 



In this chapter we discuss how different tools of quantum transport theory, based on 
the nonequilibrium Green function (NEGF) techniquesp can be extended to treat spin cur- 
rents and spin densities in realistic open paramagnetic semiconductor devices out of equi- 
librium. In such devices, the most important spin- dependent interaction is the SO coupling 
stemming from relativistic corrections to the Pauli-Schrodinger dynamics of spin-| electrons. 
Recent theoretical efforts to understand spin transport in the presence of SO couplings by 
using conventional approaches, such as the bulk conduct i vity o f infinite homogene ous sys- 



tems [computed via the Kubo formula (Murakami et al.l . j2003l ; ISinova et all 120041 ) or the 



kinetic equation ( Mishchenko et al. . 2004 )] or spin-densitj|§ diffusion equations for bounded 



both storage and information processing on a single chip. 

2 The NEGF techniques for finite-size devices have been developed over the past three decades mainly 



throu gh the studies of charge curr ents of non-interacting quasiparticles in mesoscopic semiconductor (Datta, 
1995) and nanoscopic molecular (jKoentopp et all , 2008) systems where quantum-coherent effects on elec- 
tron transport are the dominant mechanism because of the smallness of their size. Current frontiers of the 
NE GF theory are also concerned with the inclusion o f many-particle interaction effects responsible for dephas- 
ing ( Okamotol . l2007tlThygesen and and Rubiol . [2008h . For a lucid introduction to the general scope of NEGF 
formalism applied to finite-size devices see the chapter 1 by S. Datta in this Volume, as well as the chapter 
23 by K. S. Thygesen and A. Rubio in the same Volume focusing on the inclusion of electronic correlations in 
NEGF applied to molecular junctions. 

3 Spin densities within the sample and spin ac cumulation along the sample boundar i es are typically probed 
in recent SHE experiments on semiconductors (jKato et al. . 2004 d Wunderlich et al. . 2005 : Sih et all 2005 . 
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systems (IBurkov et all 12004 ; iBleibauml . 120061 ; iGalitski et all 120061 ) . have encountered enor- 
mous challenges even when treating non-interacting quasiparticles. Such intricacies in sys- 
tems with the intrinsic SO couplings, that act homogeneously throughout the sample, can be 
traced to spin non-conservation due to spin precessio n which leads to ambiguity in defining 



spin currents (jShi et all 120061 : ISugimoto et all 12006) in the bulk or ambiguity in supplyi ng 



the boundary conditions for the diffusion equations (iBleibauml . 12006c IGalitski et all 120061 ) 



On the other hand, spin-resolved NEGF techniques discussed in this chapter offer con- 
sistent description of both phase- coherent (at low temperatures) and semiclassical (at finite 
temperatures where dephasing takes place) coupled spin and charge transport in both clean 
and disordered realistic finite-size devices attached to external current and voltage probes, as 
encountered in experiments. The physical quantities that can be computed within this frame- 
work yield experimentally testable predictions for outflowing spin currents, induced voltages 
by their flow, and spin densities within the device. The presentation is tailored to be mostly 
of a tutorial style, introducing the essential theoretical formalism and practical computational 
techniques at an accessible level that should make it possible for graduate students and non- 
specialists in physics and engineering to engaged in theoretical and computational modeling 
of nanospintronic devices. We illustrate formal developments with examp l es drawn from the 



coupled semiconductor nanostructures 



filed of the mesoscopic SHE (Nikolic et all 12005a ISheng and Tina. 12006c lHankiewicz et al 



2004c iRen et all 12006c iBardarson et all 20071 ; ISilvestrov et. all 120091 ) in low-dimensional SO- 



2 WHAT IS PURE SPIN CURRENT? 

Pure spin current represents flow of spin angular momentum which is not accompanied by 
any net charge transport [Fig. (He)]. They can be contrasted with traditional electronic 
circuits where equal number of spin-f and spin- J, electrons propagate in the same direction, 
so that total charge current in that direction J = I* + 1^- is unpolarized 1=0 [Fig. QJa)]. 
Spin currents are substantially different from familiar charge currents in two key aspects: 
they are time-reversal invariant and they transport a vector quantity. In metal spintronic 
devices, ferromagnetic elements polarize electron spin thereby leading to a difference in charge 
currents of spin-f and spin- J. electrons [Fig. LUb)]. Such spin-polarized charge currents are 
accompanie d by a net spin curren t I s = -^(1^ — 1^) ^ 0, as created and detected in magnetic 
multilayers ( Maekawa and Shinjo . 20021 ). Figure [T] also provides a transparent illustration of 



one of the major advantages of pure spin currents over the spin-polarized charge currents 
employed by the "first-generation" spintronic devices. For example, to transport information 
via 2 x h/2 spin angular momenta, Fig. [TJ(b) utilizes four electrons. In Fig. [T^c) the same 
transport of 2 X h/2 is achieved using only two electrons moving in opposite direction. The 
Joule heat loss in the latter situation is only 25% of the dissipative losses in the former case. 



2006). 

4 Afi 



After learning about the NEGF techniques for spin transport discussed in this chapter, an interested reader 



articles (Souma and Nikolid 


2004; Nikolic and Soumal 2005a: Nikolic et all 2005ti ISouma and Nikolic 


2005: 


Nikolic et al.. 2005c 


. 2006: Nikolic and Zarbo. 2007: Draeomirova and Nikolic. 2007: Zarbo and Nikolic 


2007: 


Drascomirova et. al. 


2008: Nikolic and Draeomiroval 20091: IChen et al. 12009). 



6 



HOW CAN PURE SPIN CURRENTS BE GENER- 
ATED AND DETECTED? 



Among a p l ethor a of imaginative theoretical proposals (jSharmal . 120051 : iTserkovnyak et al 



20051 ; iTangj . 120061 ; iNagaosal . 120081 ) to generate pure spin currents using quantum effects in 



ferromagnet, semiconductor, and superconductor systems and their hybrids, only few have 
received con tinuous experimental attention. These include: non-lo cal spin injection in lateral 



spin valves (IValenzuela and Tinkhaml . 120061 ; iKimura et all 120071): adiabatic quantum spin 



pumps based on semiconductor quantum dots (IWatson et all 120031 ); spin current pumping 



by precessing magnetization of a ferromagne t ic lay e r driven by mic r owaves under the ferro- 



magn etic resonance conditions (ISaitoh et all 120061 ; I Cost ache et all 120061; iMorivama et. al. 



20081 ); optical pump- probe experiments on semiconductors (IStevens et al 
SHE flSih et alllioOffl . 



20031 ): and the 



Even if spin currents are induced easily, their detection can be quite challenging since trans- 
port of electron spin between two locati ons in real space is alie n to Maxwell electrodynamics 



and no "spin current ammeter" exists (lAdagideli et all 120061). In m etal spintronic devices 



spin currents can be converted into voltage signal (jJung and Led . 120051 ) by injection into ferro- 
magnetic electrode (as achieved in lateral spin valves). On the other h and, for semiconduc tor 
spintronic devices, which do not couple well to metallic ferromagnets (IFabian et all 120071 ). it 
is important to avoid ferromagnetic elements (and their stray fields) in both the spin injection 
and the spin detection processes. Multifarious theoretical ideas have been contemplated to 
solve this fundamental problem, ranging from the detection of tiny electr ic fields induced by 



the flow of magnetic dipoles associated with spins (IMeier and Losd . 20031) to nanomechanica l 
detection of oscillations induced by spin currents in suspended rods (jMaFshukov et all 120051 ) . 
Desirable schemes to detect pure spin current in semiconductors should exploit fundamental 
quantum-mechanical effects that can transform its flux into conventionally measurable voltage 
drops and charge currents within the same circuit through which the spin current is flowing. 
The recently discovered SHE holds a great promise to revolutionize generation, control, and 
detection of pure spin fluxes within the setting of all-electrical circuits. 



4 WHAT IS THE SPIN-HALL EFFECT? 



The SHE actually denotes a collection of phenomena manifesting as transverse separation of 
spin-f and spin-j states driven by longitu dinally injected standard unpqlarized charge cur 



rent o r longitudin a l external electric fiel d ([Murakami 120061 : ISchliemannl . 120061 ; ISinova et al. 



20061 ; lEngel et all 120071 ; INagaosal . 120081 ) . The spins separated in this fashion comprise ei- 
ther a pure spin current or accumulat e at the later a l sample boundaries. Its Onsager recip 
rocal phenomenon — the inverse SHE (jHirschl . Il999l ; lHankiewicz et all 12004 . 120051 ; iLi at al. 



20061 ) where longitudinal pure spin current generates transverse charge current or voltage 
between the lateral boundaries — offers one of the most efficient schemes to detect elusive 



pure spin currents by converting them into e lectrical quantities ( IValenzuela and Tinkham 



2006; 


Saitoh et al.. 


2006 





Kimura et all 120071 ). For an illustration of possible experimental 



manifestations of the direct and inverse SHE in multiterminal nanostructures see Fig. [2J 
While SHE is analogous to the classical Hall effect of charges, it occurs in the absence of any 
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Figure 2: Basic phenomenology of the direct and inverse SHEs: (a) conventional (unpolarized) 
charge current flowing longitudinally through the sample experiences transverse deflection of 
opposite spins in opposite direction due to SO coupling induced "forces" . This generates pure 
spin current in the transverse direction or spin accumulation (when transverse electrodes 
are removed) of opposite sign at the lateral sample edges; (b) pure spin current flowing 
through the same sample governed by SO interactions will induce transverse charge current 
or voltage drop AV = V2 — V3 (when transverse leads are removed). Note that to ensure purity 
(I2 = I3 = 0) of the transverse spin-Hall current in (a), employed to define manifestations of 
the mesoscopic SHE in ballistic or disordered SO-coupled multiterminal nanostructures, one 
has to apply proper voltages V 2 and V 3 . 



externally applied magnetic fields or magnetic ordering in the equilibrium state. Instead, both 
the direct and the inverse SHE essentially require the presence of some type of SO interactions 
in solids. Although SO couplings are a tiny relativistic effect for electrons in vacuum, they 
can be enhanced by several orders of magnitude for itinerant electrons in semic onductors 
due t o the interplay of crystal symmetry and strong electric fields of a tom cores flWinkler) . 



20031 ). They have recently emerged as one of the central paradigms ( IFabian et all 120071 ) 



of semiconductor spintronics — unlike cumbersome magnetic fields, they make possible spin 
control on very short length and time scales via electric fields, and could, therefore, enable 
smooth integration with conventional electro nics. 



Some of the observed SHE manifestations (IKato et all 12004 d ) have been explained (jEngel et al. 



20071 ) by the SO coupling effects localized to the region around impurities which bring inter- 



play of skew-scattering (asymmetric SO-dependent scattering which deflects spin-j and spin-], 
electrons of an unpolarized flux in opposite directions transverse to the flux) and side jump 
[due to the noncanonical nature of the p hysical position and velocity operators in the presence 
of the SO coupling around an impurity (jSinitsynl . 120081 )]. These impurity- driven mechanisms 
were crucial ingredient of the seminal arguments in early 1970s (ID'yakonov and Perell . Il97ll ) 
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predicting theoretically the existence of [in modern terminology (IHirschl . Il999l )] the extrinsic 
SHE. The extrinsic SO effects are fixed by the materials properties an d the corresponding SHE 
is ha rdly controllable, except through charge density and mobility (jAwschalom and Flatte . 



20071 ). 



A strong impetus for the reviva l of interest into th e realm of SHE has ascended from 
speculations (Murakami et all 120031 : ISinova et all 120041 ) that transverse pure spin currents, 
several orders of magnitude larger than in the case of extrinsic SHE, can be driven by longi- 
tudinal electric fields in systems with intrinsic SO couplings. Such SO couplings manifest in 
materials with bulk inversion asymmetry or semiconductor heterostructures where inversion 
symmetry is broken structurally. They act homogeneously throughout the sample inducing 
the spin-splitting of the quasiparticle energy bands. The intrinsic (or ban d structure-driven) 



SHE could account for large spin current sig nals observed in 2 D hole gases (jWunderlich et al. 



20051 ) or huge SHE response in some metal s (IGuo et al.l . l2008l ). In addition, since the strength 



of the Rashba SO coupling (jWinklerl . 120031 ) in two-dimensional electron gases (2DEGs) within 
heterostructures wi th strong structural inversion asym metry can be controlled experimentally 
by a gate electrode ( Nitta et al. . 1997 : Grundler, 2000l ). intrinsically driven SHEs are amenable 



to eas y all-electrical manipulation in realistic nanoscale multiterminal devices (jNikolic et al 



2005J). 



The magnitude of both the extrinsic and intrinsic SHE also depends on the impurities, 
charge density, geometry, and dimensionality. Such a variety of SHE manifestations poses im- 
mense challenge for attempts at a unified theoretical description of spi n transport in th e pres- 
ence of relativistic effects. This has not been resolved by early hopes (IMurakami et all 12003 ; 



Sinova et all 120041 ) that auxiliary spin current density computed within infinite homogeneous 



systems could be elevated to universally applicable and experimentally measurable quantity 
(for more technical discussion of these issues see Sec. [7]). Thus, theoretical analysis has increas- 
ingly been shifted toward experimentally relevant quantities in confined geometries and p redic- 
tions on how to control parameters that can enhance them (Onoda and Naeaosa . 2005a). Ex- 
ample s of such quantit i es are edge spin accumulat i on (Nikolic et al. . 2005 J ; Onoda and Nagaosa . 



2005a iNomura et all 120051; IZvuzin et all 120071: ISilvestroy et. al .. 2009) and bulk spin den- 



sity (Nikolic et alll2006l : iRevnoso et all 1200ft iFinkler et all 12 007: Ch en et all 120071 : iLiu et al. 



20071 1 , or outflowing spin currents driven by them (jNikolic et all 120061 ) . 



5 WHAT IS THE MESOSCOPIC SPIN-HALL EFFECT? 



Realistic devices on which SHE experiments are performed are always in contact with ex- 
ternal electrodes and circuits which typically inject the charge current (rather than applying 
"longitudinal electric field" ) or perform measurement of the resul ting voltage drops a nd spa- 
tial distribution of spins and charges. In the seminal arguments ( Sinova et all 2004 ) for the 
intrinsic SHE in infinite Rashba spin-split 2DEG (in the "clean" limit), electric- field-driven 
acceleration of electron momenta and associated pr ecession of sp i ns play s a crucial role. On 



the other hand, mesoscopic SHE was introduced (jNikolic et all 12005 a ) for ballistic finite- 
size 2DEGs attached to multiple current and voltage probes where electric field is absent 
in the SO-coupled central sa mple (on the proyiso t hat surrounding leads are reflectionless). 



Another stunning difference (jSheng and Ting . 120061 ) between intrinsically driven SHE in the 
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bulk and finite-size 2DEGs is extreme sensitivity to disorder in the former case (jlnoue et al. 



2004 ) which, for linear in momentum SO couplings (such as the Ras hba one) , is able to com- 
pletely destroy the spin-Ha ll current density in unbounded systems (IMishchenko et al.L 12004 ; 



Adagideli and Bauerl . 120051 ). Unlike in three-dimensional semiconductor and metallic devices, 



which are always disordered and where extrinsic contribution to the SHE is therefore present 
or dominant, ballistic conditions for the mesoscopic SHE can be achieved in low-dimensional 
semiconductor systems. In fact, the very recent experiment on nanoscale H-shaped structures 
realized using high mobility HgTe/HgCdTe quantum wells has reported for the first time the 



detec tion of mesoscopic SHE via non-local and purely electrical measurements (IBriine et al. 



2008|). 



The magnitude of pure spin currents flowing out of mesoscopic SHE device (illustrated 
in Fig. [5]) through ideal (spin and charge interaction free) electrodes is governed by the spin 
precession length Lso This mesoscopic length scale (e.g., L$o ~ 100 nm in recently fabri- 
cated 2DEGs), on which the vector of the expectation va lues of spin precess es by an angle tt, 



has been identified through intuitive physical arguments (lEngel et al.L 120071 ) as an important 
parameter for spin distributions (e.g., in clean systems the spin response to inhomogeneous 
field diverges at the wave vector q = 2/L so ). In fact, the mesoscopic SHE analysis pre- 



dicts ( INikolic et al.L 12005 a ) via numerically exact calculations (see Fig. |S|) that optimal device 
size for achieving large spin currents is indeed L ~ -^so- T his is further confirmed by a lter- 
native analyses of the SHE response in finite-size systems (iMoca and Marinescul . 120071 ). In 
the general cases f Sih et al. . 2005 : Hankiewicz and Vignale . 20081 ). where both the extrinsic 
and intrinsic SO interaction e ffects are present, the int rinsically driven contribution to SHE in 
finite-s ize device s dom inates ( INikolic and Zarbol . 120071 ) when the ratio of characteristic energy 
scales ( jNagaosal . 120081 ) for the disorder and SO coupling effects satisfies A so r/h > 1 (A so 
is the spin-splitting of quasiparticle energies and H/r is the disorder induced broadening of 
energy levels due to transport scattering time r). 

For mesoscopic SHE devices in the phase-coherent transport regime (device smaller than 
the dephasing length), one can also observe the effects of quantum confinement and quantum 
interferences in sp in-rela t ed qu antities that counterpart familiar examples from mesoscopic 



charge transport (jDattal. 119951 ). They include: SHE conductance fluctuations (IRen et al. 



2006; Bardarson et al. 



7l2007h: resonances in SHE conduct ance due to opening of new conduct 



ing ch annels (INikolic et al.l . l2005a : ISheng and TingJ . 120061 ) or mixing of bound ( IBulgakov et al. 



19991 ) and propagating stat es due to SO couplings; an d constructive or destructive quantum 



interference-based control (lSouma and Nikolid . 120051 ) of spin-Hall current in multiterminal 
Aharonov-Casher rings (as the electromagnetic dual of Aharonov-Bohm rings where SO cou- 
pling, rat her than magnetic field, permeates the ring). The charge and spin dephasing can be 
included ( iGolizadeh-Mojarad and Dattal . 12007a ) within the same NEGF transport formalism 
to allow for comparison with expe riments performed at finite temperatu res where quantum 
coherence effects are smeared out (IGolizadeh-Mojarad and Dattal . l2007fl ) . 

In general, the presence of SO couplings requires to treat the whole device geometry when 
studying the dynamics of transported spin densities. For example, the decay of nonequilib- 
rium spin polarizations in balli stic or disordered quantum wires is highly depend ent on the 



transverse confinement effects (jNikolic and Soumal. 12005a; iHolleitner et al.L 120061 ) or chaotic 



vs. regular boundaries of quantum dots (I Chang et all 120041 ) . Since SO couplings in SHE 
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devices manifest through both of their aspects — creation of spin currents and concurrently 
relaxation of spins — it is a nontrivial task to understand how spin c urrents and edge spin 
accum ulations scale with increasing the strength of the SO couplings (jOnoda and Nagaosa . 



2005 d ). 



The analysis of the whole device setup, where central SO-coupled sample is treated to- 
gether with the surrounding el ectrodes, also simplifies the d iscussion of esoteric SHE concepts, 
such as the SHE in insulators ( Onoda and Nagaosa . 200561 ) [where electrodes introduce dissi- 
pation necessary t o obtain nonzero y alue of spin accumulation as time-reversal odd quantity] 
or quantum SHE ( Konig et al. . 20081 ) whose quantized spin-Hall conductance is due to chiral 
spin- filtered (or helical) edge states (i.e., Kramers doublets of states forcing electrons of op- 
posite spin to flow in opposite directions along the edges of the sample) in a m ultiterminal 



SO-coupled bridge with energy gap in t he central sample (ISheng et all 120051 ) . For exam 



pie, recent direct experimental evidence (IRoth et all 120091 ) for nonlocal transport in HgTe 
quantum wells in the quantum spin-Hall regime, which shows how non-dissipative quantum 
transport occurs through chiral spin-filtered edge states while the contacts lead to equilibra- 
tion between the counter-propagating spin states at the edge, can be analyzed only by the 
device-oriented quantum transport techniques discussed in this chapter. 



6 SO COUPLINGS IN LOW-DIMENSIONAL SEMI- 
CONDUCTORS 



The coupling between the orbital and the spin degree of freedom of electrons is a relativistic 
effect described formally by the nonrelativistic expansion of the Dirac equation in external 
electric and magnetic fields (for which exact solutions do not exist) in po wers of the inverse 
speed of light c. In the second order v 2 /c 2 , one identifies (Zawadzki, 2005 ) the SO term@ 



Hi 



h 



so 



Arrive 2 



p-(ax W(r)), 



(1) 



responsible for the entanglement of the spin and orbital degrees of freedom in the two- 
component nonrelativistic Pauli Hamiltonian for spin-| electron. Here mo is the free electron 
mass, a = (a x ,a y ,a z ) is the vector of the Pauli matrices, and V(r) is the electric poten- 
tial. The SO coupling term can also be extracted from the semiclassical analysis that usually 
invokes the interaction of the electron magnetic d ipole moment (associated with spin) with 
magnetic field in the frame moving with electron (Jackson, 19981 ). In the instantaneous rest 
frame of an electron, magnetic field is obtained by Lorentz transforming electric field from the 
laboratory frame. It is actually m ore efficient for intuitive analysis of different experimental 
situations to remain (jFisherl . Il97ll ) in the laboratory frame where a magnetic dipole /j, moving 
with velocity v generates electric dipole moment 



Piab = v x n/c 2 



(2) 



5 Although a topic of numerous textbooks on relativisti c quantum mecha nics and quantum field theory, the 
v 2 /c 2 expansion has recently been carefully reexamined ( Zawadzki , 20051 ) to find all terms at this order of 
approximation in manifestly gauge invariant form, thereby revealing various inconsistencies in the textbook 
literature. 
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Here the right-hand side is evaluated in the electron rest frame and Pi a b is measured in the 
lab (both sides can be evaluated in the lab yielding the same result to first order in v/c). The 
potential energy of the interaction of the electric dipole with the external electric field Ei ab in 
the lab frame, C/di po ie = — Piab - Ei ab , corrected to the Thomas precession^] f/rhomas = — C^dipoie/2 
leads to the SO coupling term U so = U dipo \ e + U Thomas = -Pi ab • Ei ab /2. 

The lab frame analysis allows one to quickly depict the Mott skew-scattering off a target 
whose Coulomb field deflects a beam of spin-j and spin-j p articles in opposite directions, 
thereby, e.g., polarizing the beam of neutrons (Fisher . 19711) or generating skew-scattering 
contribution (jEngel et all 120071 : lHankiewicz and VignaleL 120081 ) to the extrinsic SHE. For 
example, if we look at spin-j electron from behind moving along the y-axis, whose expectation 
value of the spin vector is oriented along the positive z-axis so that the corresponding magnetic 
dipole moment lies along the negative z-axis, then in the lab frame we also see its Lorentz 
transformed electric dipole moment Pi ab oriented along the negative x-axis. The electric dipole 
feels the force F = (P lab - V)E lab , oriented in this case along the positive x-axis (right transverse 
direction with respect to the motion of the incoming electron) since gradient of the electric 
field Ei a b generated by the target is always negative outside of it. Note that this simple-minded 
classical picture only explains one aspect of the SO-dependent interaction with impurity. The 
other one — the so-called side jump (i. e., sideways sh ift of the scattering wave packet) — requires 
more quantum mechanical analysis ( Sinitsyn . 20081 ) to extract additional contribution to the 
velocity operator due to impurity potential Vdi S order( r ) i n the SO Hamiltonian Eq. (CQ). 

The heuristic discussion based on the Lorentz transformatio ns gives only a min uscule ef- 
fect and the influence of electronic band structure is essential (jEngel et all . 120071 ) to make 
these effects experimentally observable in solids. In the case of atoms, SO coupling is due to 
interaction of electron spin with the average Coulomb field of the nuclei and other electrons. 
In solids, V(r) is the sum of periodic crystalline potential and an aperiodic part contain- 
ing potentials due to impurities, confinement, boundaries and external electric fields. The 
nonrelativistic expansion of the Dirac equation can be viewed as a method of systematically 
including the effects of the neg ative-energy solu tions on the positive energy states starting 
from their nonrelativistic limit (jZawadzkil . 120051 ). This effect in vacuum is small due to huge 
gap 2rriQC 2 between positive and negative energy states. In solids, strong nuclear potential 
competes with this huge denominator in Eq. (OQ), so that much smaller band gap between 
conduction and valence band (playing the role of electron positive energy sea and positron 
negative energy sea, respectively) replaces 2mnc 2 , thereb y illustrating the origin of strong 
enhancement of the SO couplings in solids (jWinklerl . 120031 ). 

Although intrinsic SO couplings can always be written in the Zeeman form cr ■ Bso(p), 
their effective magnetic field Bso(p) is momentum-de pendent and, the refore, does not break 
the time-reversal invariance. The Kramers theorem ( jBallentind . Il998l ) for time-reversal in- 
variant quantum systems requires that the energy bands e n (k) of an electron in a periodic 
potential satisfy e n (k, j) = e n (— k, J.) since k i— ► — k and cr=ji— > er=j upon time reversal 
(Kk is crystal momentum). Therefore, in semiconductors invariant under spatial inversion 
k i— > — k (such as silicon) the Kramers theorem gives double degenerate spin states for any 
k value, e n (k, j) = e n (k, j). To obtain a non-zero Bso(p) that breaks the spin degeneracy 



6 The Thomas precession takes into account change in rotational kinetic energy due to the precession of the 
accelerated electron as seen by laboratory observer in the "extended" special relativity of accelerated objects. 
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in three-dimensional crystals £ n (k, |) 7^ £ n (k, |) the host crystal has to be inversion asym- 
metric. In bulk semiconductors with zinc-blende symmetry, the conduction band of III-V 
compounds will split into two subbands where anisotropic spin-splitting is proportional to k 3 . 
Such SO-induced splitting is termed cubic Dresselhaus SO coupling, an d it is as s ociate d with 
bulk inversion asymmetry (BIA). In semiconductor heterostructures (I Winkle rl . I2OO3L ideal 
symmetry of the 3D host crystal is broken by the interface where 2DEG or two-dimensional 
hole gas (2DHG) is confined within a quantum well. In such reduced effective dimensionality, 
the symmetry of the underlying crystal lattice is lowered, so that an additional linear in k 
Dresselhaus term becomes relevant. Besides microscopic crystalline potential V crysta \(r) as the 
source of the electric field in the SO coupling term in Eq. (TjQ), an interface electric field ac- 
companying the quantum well structural inversion asymmetry (SIA) gives rise to the Rashba 
spin-splitting of conduction band electrons 



Rashba 



a 

n iaxp) 



(3) 



for 2DEG in the xy-plane and e z as the unit vector along the z-axis. In narrow gap semi- 
conductors the Rashba effect linear in k should dominate over the bulk k 3 Dresselhaus term. 
Moreover, it has been experimentally dem onstrated that Rashba coupling ca n be changed by 



as much as 50% by external gate electrode (jNitta et al.l . 119971 : iGrundlerl . l2000l ) covering 2DEG, 
which has become on of the key concepts in semiconductor spintronics. Nevertheless, there 
is a lengthy theoretical debate on the importance of different electric field contributions to 
the value of experimentally observed a and t he parame t ers w hich are effectively manipulated 



er 



2000) —we refer to a comprehensive 



Winkler! (120031 ). Note that impurity de- 



via the gate electrode to c ause its increase fjGrund 
overview of these issues bv iFabian et all and 
termined Vdi SO rder(r) contribution to Eq. (CQ) does not require broken inversion asymmetry of 
the pure crystal or of the structure. 

The coupling between electron mom entum and sp in correlates charge currents and spin 
densities in SO-coupled semiconductors (jSilsbed . I2004T ) leading to highly non-trivial effects in 
nonequilib rium situations. Some of these have been observed i n recen t magneto-electric ex- 
periments flGanichev et al.l . l200fil : ISilov et all . I2004J : iKato et al.l . l2004fih. where charge curren t 
induces spin density, as well as in the spin galvanic experiments ( jGanichev and Prettll . 120031 ) , 
where nonequilibrium spin density drives a charge current. It is also the key ingredient of 
intrinsically driven SHEs [where the induce d spin density is oriented out-of-plane, rather than 
in-plane as in the magneto-electric effects (jSilsbed . 120041 )]. 

In this Section we focus on the description of the Rashba SO coupling^ in the form suitable 
for spin-dependent NEGF calculations using local orbital basis and the lattice Hamiltonian 
defined by it. The Rashba Hamiltonian prepared in this fashion will be used as the start- 
ing point for illustrating SHE-related spin and charge transport calculations in Sec. [HJ We 
also add the treatment of the extrinsic SO coupling in 2DEG within the same local orbital 
basis fram ework in Sec. 16.3.21 to enable the description of the most general experimental 



situations ( Sih et all 120051 ) in low -dimensional devices wh e re both extrinsic and in t rinsic 



mechanisms can ac t conc urrently (jNikolic and Zarbol . 120071 ; lHankiewicz and Vignald . 12008 



Dragomirova et. all 120081 ) 



7 The same analysis applies to linear Dresselhaus coupling since they can be transformed into each other 
by a unitary matrix. 
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6.1 Rashba coupling in bulk 2DEG 

The effective single-particle SO Hamiltonian for a clean infinite homogeneous 2DEG with the 
Rashba coupling Eq. (jHJ) can be formally rewritten as 



2D 
R 



a 



2m* 



Is + V (Py ® 0"i 



(4) 



Here <E> stands for the tensor product of two operators acting in the tensor product Ho <8> 7~Ls 
of the orbital and spin Hilbert spaces, and m* is the effective mass. To enforce pedagogical 
notation, we also use Io as the unit operator in Tio and Is for the unit operator in Kg. 
The internal momentum- dependent magnetic field corresponding to the Rashba coupling is 
extracted from Eq. (jHJ), recast in the form of the Zeeman term — g/iB<x • Br(p)/2, as Br(p) = 
(2a/gfiB) (p x e z ). The Hamiltonian commutes with the momentum operator p, time-reversal 



H,p 







operator T, and the chirality operator (<r x p/|p|) • e z . The zero commutator 

due to the translation invariance of infinite 2DEG implies that solutions of the Schrodinger 
equation are of the form 

*(r) = Ce* kr | X > s , (5) 



H,T 
E(-k 



0, the 



-o- 



The 



where \x) s is a two-component spinor. Due to time-reversal invariance 

two eigenstates of opposite momenta and spins are degenerate, -E^k, a) = 
time-reversal op erator can be writ ten as T = K ^exp(ina y /2), where K is the complex conju- 
gation operator (jBallentinei . Il998l ) . If the Hamiltonian is invariant to chirality transformation, 
the spin state is momentum-dependent. In the Rashba case the spin and momentum of an 
eigenstate are always perpendicular to each other. 

The Rashba Hamiltonian Eq. (j4j) in the case of a ID electron gas (1DEG) reduces to 



- — - - ock x a. 
2m* 



Its eigenvalues define the energy-momentum dispersion 



with the corresponding eigenvectors 



2m* 



(6) 



(7) 



e ik x x l I 1 



(8) 



These eigenstates are labeled with the momentum operator eigenvalue hk x and the chirality 
operator eigenvalue A = ±1. Equation (jHJ) shows that 



(4>k x ,±\v y \^k x ,±) = =F 



(9) 



meaning that each spin-split parabolic subband has a well-defined spin. 
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Figure 3: The one- dimensional energy-momentum dispersion for: a) infinite clean 1DEG, b) 
1DEG in magnetic field, and c) 1DEG with the Rashba SO coupling. The SO-induced spin 
splitting in (c) is signifies by the energy separation Ago at a given k x . The complementary 
description of the SO-induced spin splitting is k^o difference between the momenta of two 
electrons of opposite spins and chiralities in panel (c). 



To elucidate the physical meaning of the SO-induced spin splitting, we compare in Fig. [3] 
the Rashba dispersion with more familiar Zeeman splitting of 1DEG placed in the external 
(momentum- independent) magnetic field. In the case of the Zeeman splitting, the spin- j and 
spin- 1 subbands are shifted vertically with respect to each other by the Zeeman energy Ag. 
On the other hand, the Rashba spin-splitting depends on momentum 

Aso(^) = 2ak x , (10) 

so that E + (k x ) and E_(k x ) are shifted horizontally along the momentum axis rather than 
along the energy axis. This ensures that the system remains spin-unpolarized in equilibrium, 
as dictated by the time-reversal invariance and the fact that spin density is time-reversal odd 
quantity. The spin splitting can also be described using k$o as the difference between the 
momenta of two electrons of opposite spins and chiralities, and different experiments (such as 
Shubnikov-de Haas, Raman scattering, or spin precession) probe either Ago or kso illustrated 
in Fig. Etc). 

From Eq. @ we see that the y-component of spin of the eigenstate \ipk x ,x) is T f° r ^ = 
1 and k x negative, and j for A = 1 and k x positive. The opposite is true for the other 
branch E\=^i(k x ) Thus, the eigenvalues of spin a y are good quantum numbers since in ID 
fry commutes with the Hamiltonian. Following this argument, it is easy to understand the 
relative spin orientations of the states along each of the dispersion branches. 

As in ID case, the Rashba Hamiltonian in 2D commutes with momentum, chirality, and 
the time-reversal operator. Its eigenenergies 

fi 2 k 2 

^(k) = ^±«|k|, (11) 

are plotted in Fig. HI They are labeled by a 2D wave vector |k| = yjk 2 . + k 2 , and the chirality 
eigenvalues A = ±1 (i.e., the spin projection perpendicular to both k and growth direction 
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Figure 4: The energy-momentum dispersion E(k x , k y ) of the Rashba spin-split 2DEG. The 
states at the Fermi level lie on two concentric Fermi circles, E + (k) = Ep for the inner circle 
and -EL(k) = Ep for the outer circle, as shown in the right panel. The radius of the inner 
circle is ki and the radius of the outer one is denoted by k 2 . The Fermi momentum kp of the 
free particle is k\ < kp < k 2 . For a given momentum k, the spin of an electron is oriented 
either parallel [_E + (k) branch] or anti-parallel [E_(k) branch] to the momentum-dependent 
effective magnetic field Bft(k). The states of opposite momentum and spin on each Fermi 
circle are Kramers degenerate. 



along the z-axis). The corresponding eigenstates are 

gikr y I 

Finding the expectation value of the spin operator ha/ 2 in eigenstate Eq. ffT2"]) 

k k 



(12) 



(13) 



demonstrates (jWinklerl . 120031 ) that no common spin quantization axis can be found for all 
eigenstates of the Rashba spin-split 2DEG (e x and e y are the unit vectors within the xy-plane 
of 2DEG). 



6.2 Rashba coupling in quantum wires 

6.2.1 Energy dispersion of Rashba spin-split transverse propagating subbands 

The Rashba Hamiltonian for a quantum wire patterned within 2DEG (along the x-axis) 



R 2m* ' h 



+ - (a- x p) • z + Konf(y), 



(14) 
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Figure 5: The subband energy-momentum dispersion in the Rashba spin-split clean quantum 
wire. The parameter ksod (d is the width of the wire) defines the weak and strong SO coupling 
regimes. For large values of ksod the mixing of subbands is nontrivial giving rise to strong 
energy dependence of the wire conductance. 



describes quasi-one- dimensional electron gas (Q1DEG) of width d whose lateral confinement is 
accounted by the potential V con f(y). The motion along the confinement direction is quantized, 
such that the energy-momentum dispersion in the absence of SO coupling is split into spin 
degenerate quasi ID subbands with quadratic dispersion E±(k x ) = E n + h 2 k 2 /2m* in ID wave 
vector k x labeling an eigenstate of the subband with index n. The effective mass Rashba 
Hamiltonian for a quantum wire is translationally inv ariant along wire direct i on x, and can 
be separated into three terms H = H s ^ + H m i x + Hm ( jGovernale and Ziilickd . I2004J ) : 



Hid 



+ V coni (y), 



2m* 
hkso A A 

* a xVyi 



mr 



Ir 



2m 



-(k x + k so a y ) + 



h 2 k 2 
2 i 11 ^so 



2m* 



(15a) 
(15b) 
(15c) 



The term H s b defines the eigenenergies E n of the transverse confining potential, while Hm 
is the ID translationally invariant Rashba term for the quantum wire. Therefore, in the 
hypothetical case where the second term H m \ x is absent, the eigenenergies of the Hamiltonian 
Eq. (1T51) show only /^.-dependent splitting 



E ± (k x ) = E n + 



± a\k x 

2m* 



(16) 



When the expectation values -E mix of H mix between the eigenstates of H sh + ff 1D are of the 



order of AE r , 
the ratio 



E n+ i — E n , the mixing of the spin-split subbands becomes important. Thus, 
-Emix frkso nfi ( n 2 fr 2 \ 1 d ka d 



AE r: 



vnr 



d 



m 



'd 2 



^so 



7T 



(17) 



17 



gives a simple criterion to separate the strong and weak SO coupling regimes in QlDEGs. In 
the weak coupling regime L so d the subband mixing is negligible, as illustrated in Fig. 0(a). 
In quantum wires with strong Rashba effect, realized for Lso ^ d (or k$od > n) hybridization 
of quasi- ID subbands becomes important and Eq. (fl6l) can not capture the non-parabolic 
energy dispersions in Fig. E](b),(c). Unlike the weak coupling regime where the Hamiltonian 
eigenstates are also eigenstates of a y (i.e., their spins are polarized in-plane along the y-eatis), 
in the strong coupling r egime no common spin quan tum number can be assigned to states 



within a given subband (IGovernale and Ziilickd . 12004 ) . 



6.2.2 Spin precession in Rashba quantum wires 

The lifting of spin degeneracy of each transverse subband of a quantum wire by the Rashba 
coupling means that an electron injected at the Fermi energy Ep can have two different wave 
vectors within the wire, E + (k x i) = E^{k x2 ) = Ep. The wave vector k x \ labels the eigenstates 
of E + spin-split subband whose spinor part is f] \ [describing spin-f along the y-eatis], while 
k x2 labels the eigenstates of E_ subband associated with the eigenspinor ^. ) [describing 
spin-j along the y-axis]. Thus, when electron with spin-f along the z-axis is injected from 
a half-metallic fe rromagnet into the R ashba quantum wire, as proposed in the Datta-Das 



spin-FET device ( iDatta and Dad . Il990l ). the outgoing wave function from a wire of length L 
will be 

i/}(L) ~ e ik * lL ( • ) + eikx2L ( _i ) ■ ( 18 ) 
The probability to detect a particle with spin-f or with spin j along the z-axis is given by 

|((10)|^(L))|=4cos 2 ^^, (19a) 



and 



2 {kxi — k x2 )L 



|((0 1)|^(L))|=4sin ^ 2 ^ , (19b) 

respectively. Such modulation of current by spin precession with angle A# = Ak x L, where 
Ak x = k x2 — k x \ = 2m* a /h 2 stems from E+jkr?) — E_(k x -\) = 0, would provide the basis for 
the operation of the envisioned spin- FET device flDatta and Dasi . ll990h . It also introduces in 



a transparent fashion the concept of the spin precession length 

irk 2 



^so 



1m* a 



(20) 



along which the vector of the expectation values of spin precesses by an angle A# = it while 
propagating along the Rashba wire. 

We emphasize that the energy independence of Lgo holds for weak SO coupling and 
narrow wires. As discussed in Sec. 16.2.11 for strong SO coupling d > Lso, intersubband 
mixing becomes important, so that Eq. fT2"U]) is inapplicable. The L^o scale appears as the 
characteristic length scales for various processes in semiconductor spintronic systems. For 
example, in weakly disordered systems Lgo also plays the role of a characteristic scale for the 
D'yakonov-Perel' spin dephasing where nonequilibrium spin density decays due to randomiza- 
tion of Br(p) in each scattering event changing p and, therefore, the direction of magnetic field 
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around which spin precesses (IChang et all 120041 ; iFabian et al.l. 120071). It also sets a m esoscopic 



scale [~ 100 nm in heterostructures with la rge a (iNitta et al 



1997; 



GrundleiL 12000^1 at which 



one c an expect the largest SHE response (INikolic et al.l . 12005a 120061 ; iMoca and Marinescu 
20071 ). 



6.3 Discrete representation of effective SO Hamiltonians 

In order to perform numerical calculations on finite-size systems of arbitrary shape attached to 
external electrodes it is highly advantageous to discretize the effective SO Hamiltonian. Here 
we introduce a discretization scheme for the Rashba Hamiltonian. In addition, we also discuss 
discretization scheme for the extrinsic SO Hamiltonian which makes it pos sible to treat intrin- 
sic and extrinsic SO coupling e ffect on equal footing in nanostructures (INikolic and Zarbo . 
20071 : iDragomirova et. all l2008h . 

The grid used for discretization is collection of points m = (m x ,m y ) on a square lattice 
of constant a, where m x ,m y are integers. We use the following notation: ijj m = (m | ip) for 



the wave function evaluated at point m; and A r 



m 



.4 



m' > for the matrix element of 



an operator A. In finite difference methods, one has to evaluate the derivatives of a function 
f(x) on a grid {. . . , x m -i, x m , x m+1 , ...} = {..., (m - l)a, ma, (m + l)a, . . .}: 



dx 



m+l 



TO— 1 



2a 



(21) 



where f m = f(x m ). The second derivative can be computed as 

d f ~\ fm+l ^fm fm—1 



(22) 



For a particular discretization scheme, the matrix elements of the derivative operators can be 
expressed in the local orbital basis |m), which can be interpreted as each site hosting a single 
s-orbitals as in the tight-binding models of solids. For instance, using Eq. (1211) we get 



rn 



d 



dx 



n 



(m + 1 1 n) - (m - 1 1 n) <5„ im+ i - 8. 



n,m—l 



2a 



2a 



(23) 



where 8 is the Kronecker symbol. The same can be done for the operator d 2 / dx 2 using Eq. (J22 



rn 



dx 2 



n 



(m + 1 1 n) — 2 (m \ n) + (m — 1 1 n) 



'n,m-\-\ 



25n.m 8 



(24) 



n,m— 1 



It is worth mentioning that for a given grid, the differentiation operators matrix elements 
are dependent on the discretization scheme used. For example, if Eq. ( 123*1) is used to get 
d 2 /dx 2 = d/dx ■ d/dx it would lead to a different result than the one in Eq. (T241) . This is due 
to the fact that these matrix elements are not computed but approximated, and this depends 
on the selected discretization scheme. 



19 



The discrete one-particle operators are written in the second quantized notation as 

A = J2 ( m A n\ c ] m c n , (25) 



rn.n 



which, together with Eqs. Q23J), <MM, leads to 



dx 2 



-^ = ^2^ (Ccm+i - 4c m -i) , (26a) 

m 

= Yl 4 (4cm+i - 2c| ra c m + C^C m _!) . (26b) 



Equations fl26al) . fl26bl) can be used to discretize the effective mass Hamiltonian of a quasi- 
electron in the clean spin and charge interaction-free 2DEG 

n 2 h 2 ( d 2 d 2 \ 

fffa= _»I s = __(_ + _j»I s , (27) 

where the discretized version of the same Hamiltonian is 

Hf Te e ^ tmm 1 Cmcr^m' a . (28) 
mmV 

This is familiar tight-binding Hamiltonian on the square lattice with single s-orbital (r|m) = 
(f)(r — m) per site. Here c mcr (c mCT ) is the creation (annihilation) operator of an electron at 
site m = (m x ,m y ). The spin-independent hopping matrix element t mm ' is given by 



m 



h 2 d 2 



m') = \ -*° ifm = m ' ±ae ^ (29) 
1 otherwise . 



2m* dr 2 

where to — h 2 /(2m*a 2 ) is the orbital hopping parameter. 
6.3.1 Discretization of the Rashba SO Hamiltonian 

Following the general discretization procedure outlined above, we recast the Rashba Hamil- 
tonian in the local orbital basis representation as 

H = g mC mo -C mcr + ^ ] C mo .t mm /C m 'o-'. (30) 



mm' aa' 



While this Hamiltonian is of tight-binding type, its hopping parameters are non-trivial 2x2 
Hermitian matrices t m / m = (t mm f)t in the spin space. The on-site potential e m describes 
any static local potential, such as the electrostatic potential due to the applied voltage or 
the disorder which is usually simulated via a uniform random variable e m £ [— W/2, W/2] 
modeling short range isotropic scattering off spin-independent impurities. The generalized 
nearest neighbor hopping t^ m , = (t mm /) ffff / accounts for the Rashba coupling 

-t I s - itso&y (m = m' + e x ) , . 

-t l s + itsofrx (m = m' + e y ) ' 
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Figure 6: An illustration of the physical meaning of the orbital hopping to (a) and the Rashba 
SO hopping tso for spin-f (b) and spin-], (c) electrons on the tight-binding lattice. As shown 
in panel (a), the probability amplitude for a spin to hop between two sites without flipping 
is proportional to to- The z-spin \ J \) Z at the central site m in panels (b), (c) flips by jumping 
to any of the four neighboring sites with probability amplitude proportional to tso- 



through the SO hopping energy scale tso = oi/2a whose physical meaning is illustrated in 
Fig. El It is also useful to express the spin precession length Eq. f[2"U]) 



^so 



2t 



so 



(32) 



in terms of the lattice spacing a. 

A direct correspondence between the continuous effective Rashba Hamiltonian Eq. (1141) 
(with quadratic and isotropic energy-momentum dispersion) and its lattice version Eq. (1301 
(with tight-binding dispersion) is established by selecting the Fermi energy (e.g., Ep = 
— 3.8to) °f the injected electrons to be close to the bottom of the band E b = — 4.0to (so 
that tight-binding dispersion reduces to the quadratic one), and by using t G = h 2 /{2m*a 2 ) 
for the orbital hopping which yields the effective mass m* in the contin uum limit. F or ex- 
ample, the InGaAs/InAlAs heterostructure employed in experiments of Nitta et al. ( 19971 ) 
is characterized by the effective mass m* = 0.05mo and the width of the conduction band 
A5 = 0.9 eV, which sets to — = 112 meV for the orbital hopping parameter on a 

square lattice (with four nearest neighbors of each site) and a ~ 2.6 nm for its lattice spacing. 
Thus, the Rashba SO co upling of 2DEG formed in this heterostructure, tuned to a maximum 
value (jNitta et all 119971 ) a = 0.93 ■ 10~ n eVm by the gate electrode, corresponds to the SO 
hopping tso/to — 0.016 in the lattice Hamiltonian Eq. ( 1301) . 



6.3.2 Discretization of the extrinsic SO Hamiltonian 

The extrinsic SO coupling in disordered 2DEG is described by the Pauli Hamiltonian 



H 



Px+Py 



(33) 



, — - + Konf(y) + ^disorder (x, V) + A (<X X J>) • Wdisorder (x, y), 

2m* 

where A is the extrinsic SO coupling strength. The fourth term in Eq. (I3"3"j) can be rewritten 



as 



A (<X X p) • VVdisorder(a;, y) = -iX [(<9 y V disorder ) X - ((^disorder) d y ] a z . 



(34) 
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Equation (|23|) can then be used to find the discrete version of c^Vdisorder 
(m | (^disorder | m') = (m \d x \ m") (m" |Vd isorder | m') 

m" 

^disorder (m (^^x)^m' ,m+ae x ^disorder(m O'^x)^ 



(35) 



m',m— ae! 



2a 



The final expression (INikolic and Zarbd . 120071 ; iDragomirova et. al.l . 120081 ) of the discretized 
version of Eq. 



H 



m+7e. 

m,Q/3 ij v^j 



(36) 



introduces second neighbor hopping in the third term (in addition to usual nearest neighbor 
hopping in the second term). Here Aso is the dimensionless extrinsic SO coupling strength 
Aso = Xh/4a 2 and €ij Z stands for the Levi-Civita totally antisymmetric tensor denote the 
in-plane coordinate axes and dum my indices y, 7 t ake values ±1). For example, for the 2DEG 
in SHE experiments conducted by lSih et al.l (120051 ). the SO parameters are set to Aso — 0.005 
and tso — 0.003to for an additional Rashba term, assuming conduction bandwidth ~ 1 eV 
and using the effective mass m* = 0.074m . 



7 SPIN CURRENT OPERATOR, SPIN DENSITY, AND 
SPIN ACCUMULATION IN THE PRESENCE OF 
INTRINSIC SO COUPLINGS 

The theory of the intrinsic SHE in infinite homogeneous systems is formulated in terms of the 
spin current density which is not conserved in media with SO coupling and, therefore, does not 
have well-defined experimental measurement procedure associated with it. Moreover , these 



spin current densit ies can be non-zero even in thermodynamic equilibrium (IRashbal . 12003 



Nikolic et al.l . 120061 ) . The conservation of charge implies the continuity equation in quantum 



mechanics for the charge density p = e|\E'(r) 



dp ^ 
— + V 

dt 



0. 



(37) 



associated with a given wave function \l/(r). From here one can extract the charge current 
density j = eRe [ft(r)v$(r)] viewed as the quantum- mechanical expectation value [in the 
state \I/(r)] of the charge current density operator 



n(r)v + vn(r) 



(3* 



This operator can also be obtained heuristically from the classical charge current density 
j = en(r)v via quantization procedure where the particle density n(r) and the velocity v are 
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replaced by the corresponding operators and symmetrized to ensure that j is a Hermitian 
operator. 

In SO-coupled systems j acquires extra terms since the velocity operator ifvv = [f , H] is 
modified by the presence of SO terms in the Hamiltonian H. For example, for the Rashba 
SO Hamiltonian Eq. ([3]) the velocity operator is v = p/m* — (a/h)(a y e x — a x e y ). The spin 
density S l = |[\I/"'"(r)o"j\I/(r)] then satisfies the following continuity equation 



as* 



+v-r = r 



(39) 



dt ' ' J ^ 

In contrast to the charge continuity equation Eq. (1371) . this contains the spin current density 



J = TT^W n 



as well as a nonzero spin source term 



(40) 



(41) 



The nonzero F$ ^ term reflects non-conservation of spin in the presence of intrinsic SO cou- 
plings which act as internal momentum-dependent magnetic field forcing spin into precession. 
Thus, the plausible Hermitian operator of the spin current density 



Jk 2 - 



(42) 



is a well-defined quantity (a tensor with nine components where j\ describes transport of spin 
Si in the fc-direction, i,k = x, y, z) only when v is spin independent. 

The lack of the usual physical justification for Eq. (1421) in systems with intrinsic SO 



couplings leads to an arbitrariness in the defin ition of the spin curre nt (jShi et all 120061 ) 



Thus, different definitions lead to ambiguities ( jSugimoto et all 120061 ) in the value of the 
intrinsic spin-Hall conductivity a s u = jy/E x computed as the linear response to the applied 
longitudinal electric field E x penetrating an infinite SO-coupled (perfect) crystal. I t also yields 
qualitatively different conclus ions about the effect of impurities on SHE (jNagaosal . 120081 ). and 
does not allow us to connect (INomura et all 120051 ) directly the value of a s n to measured edge 
spin accumulation of opposite signs along opposite lateral edges. 

Under the time-reversal transformation, the mass, charge, and energy do not change sign, 
while the velocity operator and the Pauli matrices change sign, t — > — t => v — > — v and 
t — > — t cr — > — <r. Since the charge current density operator Eq. (IHH1) contains velocity, 
it changes sign under the time reversal t — » — t j — > — j and, therefore, has to vanish in 
the thermodynamic equilibrium [except in the presence of an external magnetic field which 
breaks time-reversal invariance thereby all owing; for circulating or di amagnetic charge cur- 
rents even in thermodynamic equilibrium (IBaranger and Stond . 119891 )]. On the other hand, 
the spin current density operator Eq. (|42[) is the time-reversal invariant (or even) quantity 
t -»• -t 

jl Jk- ^ ^ ne dock ran backward, spin current would continue to flow in the 
same direction. Thus, j\ can have non-zero exp ectation value s even in thermodynamic equi- 
librium. This has been explicitly demonstrated (jRashbal . 120031 ) for the case of an infinite clean 
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Rashba spin-split 2DEG where such equilibrium spin currents are polarized inside the plane. 
Therefore, the out-of-plane polarized sp in current density has been considered as a g enuin e 
nonequilibrium SHE- induced response (ISinova et all 12004 ; lHankiewicz and Vignald . 120081 ). 
However, within multiterminal finite-size Rashba coupled devic es even out-of-plane polarized 



200fih . This is essential information 



spin currents can have equilibrium nature (iNikolic et al.1 
for the development of a consistent theory for transport (nonequilibrium) spin currents where 
the contributions from background (equilibrium) currents must be eliminated. 

A plausible solution to these issues appears to be in defining a conserved quantity where 
non-conserved part Eq. ( HT1) is moved to the right-hand side of the c ontinuity equati on Eq. ( |39l) 
and incorporated in the definition of a new spin current operator (IShi et all 120061 ). However, 
this solution is of limited value since it cannot be used for realis tic devices with arbitrary 
boundary conditions due to sample edges and attached electrodes (INagaosal . 120081 ) . In Sec. M 
we discuss how to bypass these issues altogether by employing NEGF-based description of 
SHE in realistic finite-size devices where spin transport is quantified through: (?) total spin 
currents flowing out of the sample through attached leads in which they are conserved due to 
leads being ideal (i.e., spin and charge interaction free); (ii) local spin currents whose sums in 
the leads is equal to the total currents of (i); and (Hi) nonequilibrium spin density near the 
boundaries (i.e., edge spin-Hall accumulation that is part of the definition of SHE) or within 
the SO-coupled sample (which is not typically associated with SHE experimentally measured 
quantities, but represents natural ingredient of theoretical modeling). 



8 



NEGF APPROACH TO SPIN TRANSPORT IN MUL- 
TITERMINAL SO-COUPLED NANOSTRUCTURES 



The NEGF theory flLeeuwen et all , bond iRammerl 120071 ; lHaug and Jauhol 120071 ) provides a 
powerful conceptual framework, as well as computational tools, to deal with a variety of out- 
of-equilibrium situations in quantum systems. This includes steady-state transport regime 
and more general transient responses. Over the past three decades, it has become one of the 
major tools to study steady-state quantum transport of charge currents through small devices 
in phase-coherent regime where electron is described by a single wave function throughout 
the device. In this respect, NEGF offers efficient realization of Landauer's seminal i deas to 
accou nt for phase-coherent conduction in terms of the scattering matrix of the device (jDattal . 
19951 ). Furthermore, NEGF is much more complete framework making it possible to go beyond 
the paradigms of the Landauer-Biittiker scattering approach by including many-body effects 
in transport where electron-electron interactions, electronic correlations, and ele ctron-phonon 
interactions introduce dephasing lead ing to incoherent electron propagation ( jDattal . 



1995 



20051 : iThvgesen and and Rubid . 120081 ). 

The central concepts in NEGF description of quantum transport are: the sample Hamil- 
tonian H and its matrix representation H in a typically chosen local orbital basis (such as 
the discrete versions of the effective SO Hamiltonian discussed in Sec. 16.31) : the retarded 
Sieads arid the lesser S < self-energy matrices due to the interaction of the sample with the 
electrodes; self-energy matrices due to many-body interactions within the sample S int , S < ; 
and two independent Green functions — the retarded G and the lesser G K one. The retarded 
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Green function describes the density of available quantum-mechanical states, while the lesser 
one determines how electrons occupy those quantum states]^ 

In the case of phase-coherent calculations of total spin and charge currents, we only need 
to find the retarded spin-resolved self-energies due to the electrodes E£, ads (determining the 
escape rates of spin-cr electrons into the electrodes) and compute the retarded Green function 
elements connecting sites between interfaces where leads are attached to the sample. For 
convenience, we term this usage of NEGF "Landauer-Biittiker approach" in Sec. 18.11 since 
it allows to describe transport in multiterminal structures by constructing the transmission 
block of the scattering matrix in terms of the (portion of) retarded Green function matrix. 

For the computation of local quantities within the sample, we also need to find matrix 
elements of the retarded Green function between any two sites within the sample as well as 
G < . This requires to solve the Keldysh (integral or matrix) equation for G < . The same is 
true for incoherent transport where dephasing takes place in the sample requiring S; nt , 
self-energies to be computed (typically in the self-consistent fashion together with G and G < ), 
or for nonlinear transport where local charge density and the corresponding electric potential 
profile within the sample play a crucial role (jChristen and Buttikerl . Il996l ) in understanding 
the device current-voltage characteristics. We denote this procedure in Sec. 18.21 and Sec. 18.31 
as the "Landauer- Keldysh" approach since here the full NEGF theory is applied to finite- 
size (clean or disordered) devices attached to semi- infinite electrodes, i.e., to the Landauer 
setup where such electrodes simplify boundary conditions for electrons assumed to escape to 
infinity through them to be thermalized in the macroscopic reservoirs (ensuring steady-state 
transport) 



8.1 Landauer-Biittiker approach to total spin currents in multiter- 
minal nanostructures 

The experiments on quantum Hall bridges in the early 1980s were posing a challenge for the- 
oretical interp r etatio n of multiterminal transport measurements in the mesoscopic transport 
regim e (lAndol . 120031 ). By viewing the current and voltage probes on equal footing, iBiittiker 
(119861 ) has provided an elegant solution to these puzzles in the form of multiprobe formulas 



7 j GpqjVp - V q 



(43) 



They relate charge current I p — /J + 1^ in lead p to the voltages V q in all other leads attached 
to the sample via the conductance coefficients G pq . To study the spin-resolved charge currents 
Ip (a =|, i) of individual spin species |, J, we imagine that each nonmagnetic lead in Fig. [7] 



8 Both Green functions can be obtained from the contour ordere d Green functi o n def i ned for a n y two 
time values that lie a long Kadanoff-Baym-Keldysh time contour (jLeeuwen et al. . l2006t iRammeii 120071 : 
Haug and Jauhol |2007t ) 



9 The NEGF theory, which is often called "Keldysh formalism," (jLeeuwen et al. . 2006) has been initiated 
by the pioneering works of Schwinge r , Baym, Kadanoff and Keldy sh who considered infinite homog e neous 
syste ms out of equilibrium ( Rammer . 2007 : [Haug and Jauhol 2007 ). Its present applications ( Datta . fl995l 
l2005h to q uantum transport in finite-size systems attached to semi-infinite electrodes can be traced to an early 
analysis of ICaroli et al.l (|1971l ). 
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Figure 7: The four-terminal bridge for the detection of the mesoscopic SHE. The central 
region is 2DEG, where electrons are confined within a semiconductor heterostructure grown 
along the z-axis whose SIA induces the Rashba SO coupling. The four attached leads are 
clean, nonmagnetic, and without any SO coupling. The unpolarized (if = 0) charge current 
(I i 7^ 0) injected through the longitudinal leads induces spin-Hall current I% z = — -ff z in the 
transverse leads which act as the voltage probes V 2 — V 3 ^ 0, I 2 = I3 = 0. 



consists of the two leads allowing only one spin species to propagate (as realized by, e.g., half- 
metallic ferromagnetic leads). Upon replacement I p — > I p and G pq — > G pq , this viewpoint 
allows us to extract the multiterminal formulas for the spin-resolved charge currents 1° 
thereby obtaining the linear response relation for spin current 

/ /=^( / i- / p). ( 44 ) 

flowing through the lead p 

I p = — ^ ] [(^gp — G qp — G qp )V p — (G pq + G pq — G pq — G pq )V q ~^ . (45) 

The spin-Hall conductance of the four-terminal bridge sketched in Fig. [7] is then defined as 

r h ii h ij - ij 

Below, we simplify the notation by introducing the labels 

G hl G^ 4- G^ — G^ — G^ (A7n\ 

pq ^pq ' ^pq pq ^pgi \ ' a J 

pout _ ^tT I fjlt _ rjll _ fli (A7W\ 

pq pq pq pq pq~ \ J 

In fact, these coefficients have transparent physical interpretation: ^G^V P is the spin current 
flowing from the lead p with voltage V p into other leads q whose voltages are V q , while ^G p n q V q 
is the spin current flowing from the leads q 7^ p into the lead p. 
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The stan dard charge conductance coefficients (jButtikerl . Il986l ; iBaranger and Stond . Il989 ; 
Dattal . Il995l ) in the multiprobe Landauer-Biittiker formalism Eq. (143 j) are expressed in terms 
of the spin-resolved conductances as G pq = G^ pq + G^ pq + G^ q + G^ q . Their introduction in 
1980s was prompted by the need to describe linear transport properties of a single sample 
(with specific impurity arrangements and attached to specific probe configuration) by using 
measurable quantities [instead o f the b ulk conductivity which is inapplicable to mesoscopic 
conductors (IBaranger and StoneL Il989h ]. They describe total charge current flowing in and 
out of the system in response to voltages applied at its boundaries. 

Regardless of the detailed microscopic physics of transport, conductance coefficients must 
satisfy the sum rule G qp = ^2 q G pq in order to ensure the second equality in Eq. ( H3l) 



i.e., 

the charge current must be zero V q = const. I p = in equilibrium. On the other hand, the 
multiprobe spin current formula s Eq. (|45p app arently posses a nontrivial equilibrium solution 
V q = const. =>- Ip 7^ [found by iPareekl (120041 )] that would contradict the Landauer-Biittiker 
paradigm demanding usage of only measurable quantities. However, when all leads are at the 
same potential, a purely equilibrium nonzero term, ^{G^V P — G l pp V p ) = — Gj^Vp, 

becomes relevant for I p , canceling all other terms in Eq. (I4"5l) to ensure that no unphysical 
total s pin current if ^ can appear in the leads of an unbiased (V^=const. ) multiterminal 



device ( Souma and Nikolic . 2005 ; Kiselev and Kim . 2005 ; Scheid et al. . 20071 ) 



At zero temperature, the spin-resolved conductance coefficients 



-y 
h ^ 



t pq 



(48) 



where summation is over the conducting channels in the leads, are obtained from the Landauer- 
type formula as the probability for spin-cr' electron incident in lead q to be transmitted to 
lead p as spin-cr electron. The quantum-mechanical probability amplitude for this processes 
is given by the matrix elements of the transmission matrix t pg , which is determined onl y 
by the wave functions (or Green functions) at the Fermi energy (iBaranger and Stond . 1 19891 ) . 
The stationary states of the structure 2DEG + two leads supporting one or two conducting 
channels can be found exactly by matching the wave functions in the leads to the eigenstates 
of the Hamiltonian Eq. (fT4"j) . the reby allowing one to obtain t he charge conductance from the 
Landauer transmission formula (IGovernale and Ziilickd . 120041 ) . However, modeling of the full 
bridge geometry with two extra leads attached in the transverse direction, the existence of 
many open transverse propagating modes ("conducting channels"), and possibly strong SO 
coupling regime when sample is bigger than the spin precession length (L > L$q), is handled 
much more efficiently through the NEGF formalism] 10 ! 

For non-interacting particle which propagates through a finite-size sample of arbitrary 
shape, the transmission matrices 



t pq 

Im E„ 



y'-Im S p <g> I s ■ Gpq ■ ^-ImS 5 



2i y 1 



H) 



(49) 



10 The multiterminal an d multichannel finite-siz e device can be also modeled using the random matrix theory 
for its scattering matrix dBardarson et all l2007t) . Howe ver, this method is strictly applicable only for weak 
SO coupling regime L L§o (jAleiner and Fal'kd . 120011 ). 
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between different leads can be evaluated in a numerically exact fashion using the realcg>spin- 
space retarded Green operator G r defined in the Hilbert space Ho <8> Hs- Their matrix 
representation in a basis |m) <g> \a) G Ho <%> Hs introduced by the Hamiltonian Eq. (I3"U]) is 
obtained through the matrix inversion 

G = , (50) 

Mo®I S -H-5^= 1 S p (8)Is 

Here \a) are the eigenstates of the spin operator for the chosen spin quantization axis. The 
matrix elements (m', a' |G r |m, a) of G yield the probability amplitude for an electron to prop- 
agate between two arbitrary locations m and m' (with or without flipping its spin a during 
the motion) inside an open conductor in the absence of inelastic processes. Its submatrix G pq , 
which is required in Eq. ( )49l) . consists of those matrix elements which connects the layer of 
the sample attached to the lead q to the layer of the sample attached to the lead p. The sum 
of the self-energy matrices Ylp=i ® Is accoun ts for the "in teraction" of an open system 



with the attached four ideal semi-infinite leads p (jDattal . Il995l ) 



8.1.1 General expression for spin-Hall conductance 

Since the total charge current I p — /J + 1^ depends only on the voltage difference between the 
leads in Fig. [7J we set one of them to zero (e.g., V4 = is chosen as the reference potential) 
and apply voltage V\ to the structure. Imposing the requirement I2 = h = for the voltage 
probes 2 and 3 allows us to get the voltages V2/V1 and V3/V1 by inverting the multiprobe 
charge current formulas Eq. (I43p . Finally, by solving Eq. (|45p for Ijj! we obtain the most 
general expression for the spin-Hall conductance defined by Eq. 



r - h 



( /^yout I /^yout , /^yout\ ^2 /-i'vn ^3 ^in 
{isr 12 + Lt 32 + (J42 )— - ^23 _ Lr 21 



(51) 



This quantity is measured in the units of the spin conductance quantum e/47r (as the largest 
possible G s h when transverse leads support only one open conducting channel) , which is the 
counterpart! 11 ! of the familiar charge conductance quantum e 2 /h. 

In contrast to charge current, which is a scalar quantity, spin current has three components 
because of the vector nature of spin (i.e., different "directions" of spin correspond to different 
quantum mechanical superpositions of ||) and |j) states). Therefore, we can expect that, in 
general, the detection of spin transported through the transverse leads of mesoscopic devices 
will find its expectation values to be nonzero for all three axes. However, their flow properties 



7|« = -1^, (52a) 
Jf« = -II*, (52b) 

I S 2 V = (52c) 

show that only the z and the x components can represent the SHE response for the Rashba SO- 
coupled four-terminal bridge. That is, if we connect the transverse leads 2 and 3 to each other 



Note that (H/2e)(e 2 /h) = e/4?T, where (e 2 //i) is the natural unit for spin-resolved charge conductance 



coefficients Gp^ 
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(thereby connecting the lateral edges of 2DEG by a wire), only the spin current carrying z- and 
x-polarized spins will flow through them, as expected from the general S HE phenomenology 



where nonequilibrium spin-Hall accumulation detected in experiments (IKato et all l2004d ) 
has opposite sign on the lateral edges of 2DEG. 

Therefore, to quantify all nonzero components of the vector of tra nsverse spin current i n 



the linear response regime, we can introduce three spin conductances ( iNikolic et all 12005 a ) 



G x sH = I^/V,, (53a) 
G% = (53b) 



G z sH = It/Vu (53c) 

(assuming V4 = 0). They can be evaluated using the same general formula Eq. (l5*Tj) where the 
spin quantization axis for f, j in spin-resolved charge conductance coefficients is chosen to be 
the x, y, or z axis, respectively. For example, selecting b z \\) = + | j) and cr z \l) = for 
the basis in which the Green operator Eq. (1501) is represented allows one to compute the z- 
component of the spin current I p z . In accord with their origin revealed by Eq. (1521) . we denote 
G z sH and G x sH as the spin-Hall conductances, while G y sp is labeled as the "spin polarization" 
conductance since it stems from the pol arization of 2D EG by the flow of unpolarized charge 



current in the presence of SO couplings (jSilsbed . [200J). 



8.1.2 Symmetry properties of spin conductances 

Symmetry properties of the conductance coefficients with respect to the reversal of a bias volt- 
age or the direction of an external magnetic field play an essential role in our u nderstanding of 
linear respon se electron transport in macroscopic and mesoscopic conductors flButtikeii[l98~6l : 



Dattal . ll995h . For example, in the absence of magnetic field they satisfy G pq = G qp [which 



can be proved assuming a particular model for charge transport (jDattal . Il995l )]. Moreover, 
since the effective magnetic field B^(p) of the Rashba SO coupling depends on momentum, 
it does not break the time-reversal invariance thereby imposing the following property of the 
spin-resolved conductance coefficients G™ = G qp ~ a in multiterminal SO-coupled bridges. 

In addition, the ballistic four-terminal bridge in Fig. [7] with no impurities possesses various 
geometrical symmetries. It is invariant under rotations and reflections that interchange the 
leads, such as: {1) rotation C4 (C2) by an angle n/2 (tt) around the z axis for a square 
(rectangular) 2DEG central region; (it) reflection a vx in the xz plane; and (Hi) reflection a vy 
in the yz plane. These geometrical symmetries, together with G pq = G qp property, specify 
V2/V1 = V3/V1 = 0.5 solution for the voltages of the transverse leads when I 2 = I3 = 
condition is imposed on their charge currents. 

The device Hamiltonian containing the Rashba SO term commutes with the unitary trans- 
formations which represent these symmetry operations in the Hilbert space Tio <S> 7is ; (i) 
U(C 2 ) (g> exp (i| <j 2 ), which performs the transformation a x — > — a x , a y — > — a y , o z — > a z and 
interchanges the leads 1 and 4 as well as the leads 2 and 3; (ii) U(a vx ) ® exp (i^a y ), which 
transforms the Pauli matrices a x — > —cr x , b y — > cr y , d z — > —b z and interchanges leads 2 and 3; 
and (iii) U(a vy ) ®exp (i^a x ) which transforms a x — > cr x , a y — > —o y , b z — > — cr z and exchanges 
lead 1 with lead 4. The Hamiltonian also commutes with the time-reversal operator. 
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The effect of these symmetries on the spin-resolved charge conductance coefficients, and 
the corresponding spin conductances G z sH) G x sH and G y sp expressed in terms of them through 
Eq. (I5ip . is as follows. The change in the sign of the spin operator means that spin-f be- 
comes spin-i so that, e.g., will be transformed into — G™. Also, the time-reversal implies 
changing the signs of all spin operators and all momenta so that G^ = G~°'~° is equivalent 
to Gpg = — Gqp*. Thus, invariance with respect to U(a vy ) ® exp(ij& x ) yields the identities 
G 2 n { x = G 1 ^, G^x* = —G l 2 ^, and G l 2 { z = —G l 2 £. These symmetries do not imply cancella- 
tion of G l ^f ■ However, invariance with respect to U(a vx ) ®exp {i\d- y ) and U(C 2 ) ®exp (i| a z ) 
implies that G^ y = and G^ z = 0. 

These symmetry imposed conditions simplify the general formula Eq. ( [51]) for spin con- 
ductances of a perfectly clean Rashba SO-coupled four-terminal bridge to 

Gin = 7£{2Gg* + Gg*), (54a) 



2e 
h 



G% = fG°^\ (54b) 



G z sH = -GT Z - (54c) 

where we employ the result V2/V1 = V3/V1 = 0.5 valid for a geometrically symmetric clean 
bridge. Because this solution for the transverse terminal voltages is violated in disordered 
bridges, its sample specific (for given impurity configuration) spin conductance cannot be 
computed fro m simplified formulas Eq. as it is often assumed in the recent mesoscopic 
SHE studies flSheng and TingL 120061 : iR^net all l200^ >. 

It insightful to apply the same symmetry analysis to the bridges with other types of SO 
couplings. For example, if the Rashba term in the Hamil tonian Eq. (fl¥l) is replaced by the 
linear Dresselhaus SO term | (px& x — p y a y ) due to BIA ( Zutic et all 120041 ). no qualitative 
change in our analysis ensues since the two SO couplings can be transformed into each other 
by a unitary matrix (a x + a y )/y/2. In this case, the spin-Hall response is signified by I^ z = 
—I^ z and = —1^ components of the transverse spin current, while if* = I^ x . For the 



Dresselhaus SO-coupled bridge, the general expression Eq. ( 1511 simplifies to 

G x sp = 7^G°f' x , (55a) 
G v sH = |(2G^ + G°H, (55b) 

G z sH = Ye ™'*- ( 55c ) 

The qualitatively different situation emerges when both the Rashba and the linear Dresselhaus 
SO couplings become relevant in the central region of the bridge since in this case it is 
impossible to find spin rotation which, combined with the spatial symmetry, would keep the 
Hamiltonian invariant while only transforming the signs of its spin matrices. Moreover, for 
such ballistic bridge the condition J 2 = J 3 = leads to V 2 /Vi = 1 — V3/V1 solution for 
the voltages, whereas imposing the alternative condition V 2 — V3 generates non-zero charge 
currents flowing through the transverse leads 2 and 3 together with the spin currents [no 
simple relations akin to Eq. ( 1521) can be written in either of these cases]. 
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Figure 8: The spin-Hall conductance G z u of: (a) Rashba SO-coupled square-shaped 2DEG as 
the function its size L (in the units of the lattice spacing a ~ 3 nm) for t$o = O.Olto setting 
the spin precession length L$o « 157a; and (b) single channel Aharonov-Casher ring as the 
function of the dimensionless Rashba SO coupling Qr = (tso/to)N/ii where the number of 
lattice sites discretizing the ring is iV = 100. 



8.1.3 Example: Transverse total spin currents in mesoscopic and quantum- 
interference-driven SHE nanostructures 

Figure M plots two examples of G Z H computed for simply-connected (square-shaped) 2DEG 
and multiply-connected ring realized within 2DEG. In both cases, the Rashba SO coupling is 
present within the device (gray area in the insets of Fig. [S]). Figure [H] confirms that optimal 
device si ze to observe large m esoscopic SHE in multiterminal structures is indeed governed 
by Lso ( Nikolic et al. . 2005 th . The computation of SHE conductance for very large square 



lattices in Fig. E^a) is made possible by the usage of the recursive Green function algorithm 
for four-terminal nanostructures discussed in Sec. 19.21 

The spin conductance of the single- channel ring attached to four single - channel ideal leads 



i n Fig . [Mb) illustrates the quantum-interference driven SHE introduced by lSouma and Nikolic 



(120051 ). The ring represent a solid-state realization of the two-slit experiment since an electron 
entering the ring can propagate in two possible directions — clockwise and counterclockwise. 
The superpositions of corresponding quantum states are sensitive to the acquired topologi- 
cal phases in a magnetic [Aharonov-Bohm effect] or an electric [Aharonov-Casher effect for 
particles with spin] external field whose changing generates an oscillatory pattern of the ring 
conductance. Thus, unlike commonly discussed extrinsic and intrinsic SHE, whos e essence 



can b e understood using semiclassical arguments and wave packet propagation (ISinitsyn 



20081 ) . here the Aharonov-Casher phase difference acquired by opposite spin states during 
their cyclic evolution around the ring plays a crucial role. The spin conductance G z sH be- 
comes zero at specific values of the Rashba coupling when the destructiv e interference of op 



posit e spins traveling in opposit e direc tions around the ring takes place ( ISouma and Nikolic 



20051 : iTserkovnyak and Brataad . 120071 ). The amplitude of such quasiperiodic oscillations of 
G Z H , which are absent in simply connected mesoscopic SHE devices of Fig. El(a), gradually 
decreases at large Rashba coupling because of the reflection at the ring-lead interface. 
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8.2 Landauer-Keldysh approach to local spin currents in multiter- 
minal nanostructures 

The theory of imaging of charge flow on nanoscale can be constructed efficiently within the 
frame work of lattice mode l s of mesoscopic devi ces and the corresponding bond charge cur- 
rents (IBaranger and Stonel . Il989l ; iTodorovl . 120021 ) . This makes it possi ble to obtain a detailed 
picture of charge p r opagation between two arbi t rary s i tes of the lattice (|Nonovama and Qguri , 



19981: ICresti et all 1200.4 iMetalidis and Brunei 1200.4 INikolic et all . l2004 IZarbo and Nikolic 



2007). thereby pr oviding a way to inter pret recent scanning probe experiments. These experi 
mental advances (ITopinka et al.l . 120031 ) have brought new insights into quantum transport by 
imaging its local features within a single sample, rather than performing conventional mea- 
surement of macroscopically averaged quantities. At the same time, scanning probe techniques 
are becoming increasingly important in the quest for smaller electronic devices. For exam- 
ple, recent imaging of charge flow in conventional p — n junctions suggests that in structures 
shrunk below 50 nm individual positions of scarce dopants will affect their fu nction, thereby 



requi ring to know precisely how charge carriers propagate on the nanoscale (lYoshida et al 



20071 ). In the case of spin transport, Kerr rotation microscopy has made possible imaging of 
steady-stat e spin density driven by charge flow t hroug h various SO-coupled semiconductor 
structures (ICrooker and Smith! l2004 iKato et all 120051 ). 

In this Section, we discuss NEGF-based tools that allow us to compute the sp atial details of 



spin flow on the scale of few nanometers by introducing the bond spin currents (INikolic et al. 



20061 ). They represent the analog of bond charge currents, as well as a lattice version of the 
spin curr ent density conventiona l ly employed in the s tudies of the intrinsic SHE in macroscopic 
systems (Murakami et al.l . 120031 : ISinova et all [200J). Even though their sums over the cross 
sections within the sample change as we move from the bottom to the top transverse elec- 
trode due to spin current non-conservation, they illustrate propagation of preces sing spins and 



withi n the ideal transverse leads their sums over the cross section reproduce (INikolic et al 



20061 ) conserved total spin currents discussed in Sec. 18.11 using the Landauer-Biittiker ap- 
proach. 



8.2.1 Bond charge current operator in SO-coupled systems 

The charge conservation expressed through the familiar continuity equation Eq. (1371) yields a 
uniquely determined bond charge current operator for quantum systems described on a lattice 
by a tight-binding- type of Hamiltonian Eq. (|3"U1) . That is, the Heisenberg equation of motion 



dN n 
~dt 
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ih 



N m ,H 



for the electron number operator N m on site m, N„ 
continuity equation on the lattice 
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(56) 

, leads to the charge 



(57) 



k=x,y 



This equation introduces the bond charge-current operator (ITodorovl . 120021 ) J mm ' which de- 
scribes the particle current from site m to its nearest neighbor site m'. The 'bond' terminology 
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is supported by a picture where current between two sites is represented by a bundle of flow 
lines bunched together along a line joining the two sites. 

Thus, the spin-dependent Hamiltonian Eq. (IHUj) containing 2x2 hopping matrix defines 
the bond charge- current operator J mm ' = ^mm' which can be viewed as the sum of four 
different spin-resolved bond charge- current operators 



J GO 

mm ' ~ in 



r> t a a r — H r 
L m'cr'' / m'm c m(T 



(5* 



where H.c. stands for the Hermitian conjugate of the first term. In particular, for the case of 
^mm' being determined by the Rashba SO interaction Eq. fl3T|) . we can decompose the bond 
charge current operator into two terms, J mm ' = J mm / + ^mm'; having transparent physical 
interpretation. The first term 
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kin 
mm' 



H.c. 



(59) 



can be denoted as "kinetic" since it originates only from the kinetic energy to and does not 
depend on the SO coupling energy £ so . On the other hand, the second term 
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(60) 



represents additional contribution to the intersite charge current flow due to non-zero Rashba 
SO hopping £ so . Here we also introduce the "bond spin-density" operator 

Smm' = J + HX - ' ( 61 ) 

a/3 

defined for the bond connecting the sites m and m', which reduces to the usual definition of 
the local spin density operator for m = m' [see Eq. (1741) ]. 



8.2.2 Nonequilibrium bond charge current in SO-coupled systems 

The formalism of bond charge current makes it possible to compute physically measur- 
able ( ITopinka et all 120031 ) spatial profiles of local charge current density within the sample as 
the q uantum-statistical average^! of the bond charge- current operator in the nonequilibrium 
state ( Caroli et al. . 1971 ; iNonoyama and Oguri . 1998 ; Cresti et al. . 2003 ). 
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(62) 
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12 The quantum statistical average (A) = Tr[pA] is taken with respect to the density matrix p that has 
evolved over sufficiently long time, so that nonequilibrium state and all relevant interactions are fully estab- 
lished. 
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Here the local charge current is expressed in terms of the nonequilibrium lesser Green function 

The usage of the second quantized notation in Eq. f)3U|) facilitates the introduction of NEGF 
expre ssions for the nonequilibrium expectation values (jCaroli et al.l . ll97ll : lNonoyama and Oguri 
19981 ) . We imagine that at time t' = — oo the sample and the leads are not connected, while 
the left and the right longitudinal lead of a four-probe device are in their own thermal equi- 
librium with the chemical potentials (jll and /i^, respectively, where /i^ = fiR + eV. The 
adiabatic switching of the hopping parameter connec ting the leads and the sample generates 
time evolution of the density matrix of the structure ( Caroli et all 1971 ). The physical quan- 
tities ar e obtained as th e nonequilibrium statistical average (...) [with respect to the density 
matrix (iKeldyshl . Il965l ) at time t' = 0] of the corresponding quantum-mechanical operators 
expressed in terms of c mcr and c mcr . Th is will lead to the expressions of the type (cj^c mcr '), 



which define the lesser Green function (jCaroli et all Il97ll ; iNonoyama and Oguril . Il998l ) 



'a') 



~Gm>m,a'a( T ~ 0) 
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(64) 



Here we utilize the fact that the two-time correlation function [c mcr (t) = e lHt l h c ma e iHt / h j 



^mm'^'lM') = r (c m , CT ,(t')Cm<j(f) 



it 



(65) 



depends only on r = t — t' in stationary situations, so the time difference r can be Fourier 
transformed to energy 
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2nh 



dEG< m ,,(E)e 



iEr/h 



(66) 



which will be utilized for steady-state transport studied here. We use the notation where 
G< m , is a 2 x 2 matrix in the spin space whose a a' element is G mm , aa ,. 

The spin-resolved bond charge current in Eq. fl6"3"j) describes the flow of charges which 
start as spin a electrons at the site m and end up as a spin a' electrons at the site m' where 
possible spin-flips a ^ a 1 (instantaneous or due to precession) are caused by spin-dependent 
interactions. The decomposition of the bond charge-current operator into the kinetic and SO 
terms leads to a Green function expression for the corresponding nonequilibrium bond charge 
currents f 

J mm ' ) = ( Jmm' ) + ( ^mm' ) 5 (67) 



(68) 



with kinetic and SO terms given by 
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Tr s {[(m' - m) x <xl [G< m ,(E) + G<, m ( J E)]} . (69) 



Note, however, that "kinetic" term is also influenced by the SO coupling through G < . In 
the absence of the SO coupling, Eq. (Ifill vanishes and the bond charge current reduces to 
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the standard expression (jCaroli et all Il97ll ; iNonoyama and Oguril . Il998l ; ICresti et all 120031 ). 
The trace Trs is performed in the spin Hilbert space. Similarly, we can also obtain the 
nonequilibrium local charge density in terms of G < 
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(70) 



which is the statistical average value of the corresponding operator. 



8.2.3 Bond spin current operator in SO-coupled systems 

To mimic the plausible definition of the spin-current density operator j l k in Eq. (l4*2l) , we intro- 
duce the bond spin-current operator for the spin-S^ component as the symmetrized product 
of the spin-i operator hfri/2 (i = x, y, z) and the bond charge-current operator from Eq. (1571) 

^mm' = ^-Yl [cL' j 9{^) t m'm} /3a C ma -H.C. . (71) 

a/3 

By inserting the hopping matrix t m / m Eq. (I3T1) of the lattice SO Hamiltonian into this ex- 
pression we obtain its explicit form for the case of the Rashba coupled system 



^mm< = % \ ( £ U (*<)/*, a ma - H.c.) + t so N mm , [ e< x (m' - m)], , (72) 



a/3 

which can be considered as the lattice version of Eq. (j42l) . Here we simplify the notation by 
using the "bond electron-number operator" N mm > = \ J2 a \ $m><jC™° + H.c. j , which reduces 
to the standard electron-number operator for m = m' '. 



8.2.4 Nonequilibrium bond spin currents in SO-coupled systems 

Similarly to the case of the nonequilibrium bond charge current in Sec. 18.2.21 the nonequi- 
librium statistical average of the bond spin-current operator Eq. ( 1721) can be expressed using 
the lesser Green function G < as 

(J**) = + (73a) 

/ \ t r°° dv 

{J S m ^ n) ) = f ^^[^(G^-G^)], (73b) 
(S 0) ) = h x (m' - m)] z tf |~ ^Tr s [G< m ,(E) + G<, m (£)] . (73c) 

Here we also encounter two terms which can be interpreted as the kinetic and the SO contri- 
bution to the bond spin current crossing from site m to site m'. However, we emphasize that 
such SO contribution to the spin-^ bond current is identically equal to zero, which simplifies 
the expression for this component to Eq. ( 173bl) as the primary spin current response in the 
SHE. 
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8.3 Landauer-Keldysh approach to local spin densities 

Motivated by recent advances in Kerr rotation microscopy, which have made possible ex- 
perimental imaging of steady-state spin polarization i n various SO-coupled semiconductor 
structures (iCrooker and Smith! . 120051 : iKato et all 120051 ). we discuss in this Section a NEGF- 

based approach to computation of the spatial profiles of ^S m ^. Unlike local spin current 

density, such local flowing spin density is a well-defined and measurable quantity that offers 

insight into the spin flow in the nonequilibrium steady transport state. For ^S m ^ computed 

at the lateral edges of the sample, we use the t erm "spin accumulation" which was directl y 
measured in the seminal spin-Hall experiments (IKato et all 12004 d ; IWunderlich et all 120051 ) . 



8.3.1 Local spin density and its continuity equation 

The local spin density in the lattice models is determined by the local spin operator S r 
(5*, S^, S^J at site m defined by 
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The Heisenberg equation of motion for each component S* 
operator 



dt 
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can be written in the form 
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where J^ m , is the bond spin-current operator given by Eq. (JTTT) so that the second term on 
the left-hand side of Eq. ( 1761) corresponds to the "divergence" of the bond spin current on site 
m. Here, in analogy with Eq. (141]) . we also find the lattice version of the spin source operator 
F£* whose explicit form is 
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(77a) 
(77b) 
(77c) 



The presence of the non-zero term F~j* on the right-hand side of the spin continuity Eq. (IToT) 
signifies, within the framework of bond spin currents, the fact that spin is not conserved 
in SO-coupled systems. The fact that the bond spin current operator Eq. (ITT|) appears in 
the spin continuity equation Eq. ( 1761) as its divergence implies that its definition in Eq. ( I7T1) 
is plausible. However, the presence of the spin source operator F^ reminds us that such 
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definition cannot be made unique (IShi et all 120061 ). in sharp contrast to bond charge current 



which is is uniquely determined by the charge continuity Eq. (157)) . 

Evaluating the statistical average of Eq. (I76p in a steady state (which can be either equi- 
librium or nonequilibrium) , leads to the identity 



^m,m+e fc / \^m-e fc ,m}) ( -^m ) ■ 



k=x,y 

In particular, for the spin-S'z component we get 

((^m Z ,m+e fc ^ — (^m-e fe ,m^)) = ~ ^ ((^m,m+e fc ^) + (^m-e fe ,m^) 7 (79) 

k=x,y k=x,y 

which relates the divergence of the spin-S'z current (left-hand side) to the spin- source (right- 
hand side) determined by the sum of the longitudinal component of the spin-S^ current and 
the transverse component of the spin-S^ current. 

Since no experiment has been proposed to measure local spin current density within the 
SO coupled sample, defined through Eq. (HOI or its lattice equivalent Eq. ( 1731) . we can obtain 
additional information about the spin fluxes within the sample by computing 
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as the nonequilibrium spin density driven by charge transport in the presence of SO couplings. 

8.4 Spin-resolved NEGFs for finite-size multiterminal devices 

The spin-dependent NEGF formalism discussed in Sec. 18.21 does not actually depend on the 
details of the external driving force which brings the system into a nonequilibrium state. That 
is, the system can be driven by either the homogeneous electric field applied to an infinite 
homogeneous 2DEG or the voltage (i.e., electrochemical potential) difference between the 
electrodes attached to a finite-size sample. For example, in the latter case, the external bias 
voltage only shifts the relative chemical potentials of the reservoirs into which the longitudi- 
nal leads (employed to simplify the boundary conditions) eventually terminate, so that the 
electrons do not feel any electric field in the course of ballistic propagation through clean 
2DEG central region. The information about these different situations is encoded into the 
lesser Green function G < . 

Here we focus on experimentally relevant spin-Hall devices where finite-size central region 
(C), defined on the L x L lattice, is attached to four external semi- infinite leads of the same 
width L. The leads at infinity terminate into the reservoirs where electrons are brought into 
thermal equilibrium, characterized by the Fermi-Dirac distribution function f(E — eV p ), to 



ensure the steady-state transpor t . In such multiterminal Landauer setup (jButtikerl . 11986 



Baranger and Stond . Il989l : iDattal . Il995l ) , current is limited by quantum transmission through 



a potential profile while power is dissipated non-locally in the reservoirs. The voltage in each 
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lead of the four-terminal spin-Hall bridge is V p (p = 1, . . . , 4), so that the on-site potential e m 
within the leads has to be shifted by eV p . 

The spin-dependent lesser Green function G < defined in Eq. (1651) is evaluated within the 
finite-size sample region as a 2L 2 x 2L 2 ma trix in the site®spin space through the spin-resolved 
matrix Keldysh equation flKeldvshl . fl965h 

G < (E) = G(E)S < (E)G t (E), (81) 

which is valid for steady-state transport (i.e., when transients have died away). Within the 
effective single-particle picture, the retarded Green function is computed by inverting the 
Hamiltonian of an open system sample+leads 
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are exactly computable in the non-interacting electron approximation and without any in- 
elastic processes taking place within the sample. They account for the "interaction" of the 
SO-coupled sample with the attached leads, thereby generating a finite lifetime that electron 
spends within the 2DEG before escaping through the leads toward the macroscopic thermal- 
izing reservoirs. Here Ic is the 2L 2 x 2L 2 identity matrix and I p is the identity matrix in the 
infinite site®spin space of the lead p. We use the following Hamiltonian matrices 

(86a) 
(86b) 
(86c) 

where |l mcr ) is a vector in the Fock space (meaning that the occupation number is one for the 
single particle state \ma) and zero otherwise) and H is the Hamiltonian given in Eq. fl30l) . 

In the general case of arbitrary applied bias voltage, the gauge invariance of measurable 
quantities (such as the current-voltage characteristics) with respect to the shift of electric 
potential everywhere by a constant V c , eV p eV p + eV c and eU m eU m + eV c , is satisfied 
on the proviso that the retarded self-energies S P (_E — eV p ) introduced by each lead depend 
explicitly on the applied voltages at the sample boundary, while the computation of the 
retarde d Green function G(E) has to include the electric potential landscape^ U m within the 
sample (IChristen and Biittikerl . Il996l ). However, when the applied bias is low, so that linear 
response zero-temperature quantum transport takes place through the s ample [as determined 
by G(Ep)], the exact profile of the internal potential becomes irrelevant (jBaranger and Stond . 
19891 : iNikolic and Allenl . Il99flh . 



13 The potential landscape U m can be obtained from the Poisson equation with charge density Eq. ([TU)) as 
the source. 
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8.4.1 How to introduce dephasing into SHE nanostructures 



Much of the applications of NEGF techniques to mesoscopic semiconductor and nanoscopic 
molecular systems tre at phase- coherent t ransport by relying on some type of the Landauer 



transmission formula (jRocha et all 120061 ). On the other hand, applications demand devices 



operating at room temperature where device size is typically much bigger than the dephasing 
length so that one seldom observes quantum interference effects in the transmission due 
to multiple scattering off impurities and boundaries. Since such quantum coherence effects 
are also seen in the linear res ponse spin density in SO-coupled quantu m wires, typically 
computed at zero-temperature (jNikolic et all 12005 d ; iReynoso et all 120061 ) , to compare such 
calculations with experiments conducted at finite temperature it is desirable to introduce 
dephasing into NEGF formalism in a simple yet controllable fashion (which leaves the charge 
current conservation intact). 

The most widely used exam ple of the simp le approach to inclusion of dephasing are fic- 
titious Biittiker voltage probes (jButtikerl . Il986l ). such as one-dimensional electrodes attached 
at each site m of the lattice and with their electrochemical potential adjusted to ensure that 
the current drawn by each probe is zero0 In NEGF formalism this introduces additional 
self-energy whose matrix elements are (Sdeph)mm,o-(7 = — if} for all sites m within the sample 



except for those where the external electrodes are attached (jGolizadeh-Mojarad and Dattal . 
2007a) • However, while Biittiker probes coupled to each site can remove sharp features in the 
transmission of the device, they also introduce additional scattering events thereby artificially 
enhancing the resistance. 

A "momentum conserving" phenomenological model for deph asing, that is easily imple- 
mented within the NEGF formalism, has recently been proposed by lGolizadeh-Mojarad and Datta 
(2007a). Here the self-energies are taken as Sd ep h (E) = diG(E) and £< ph (£) = ^G<(£), 
with parameter di controlling the strength of the elastic dephasing processes. This choice 
is motivated by the general expressions for the first-order self- consistent Born approxima- 



tion for the self-energies which ensure the current conservation (ILeeuwen et all 120061 ). The 



additional choice of the function multiplying matrix elements of NEGFs is taken to guar- 
antee momentum conservation, so th a t deph asing does not artifactually enhance the resis- 
tance (jGolizadeh-Mojarad and Dattal . l2007d ). This form of the retarded and lesser self- 
energies requires that Eq. ( 182|) be solved self-consistently [with the initial guess for G(E) 
being the solution with no dephasing] and then use such G(E) to solve the Keldysh equation 
Eq. (1811) directly as the Sylvester equation of matrix alge bra. This allows one to compare 
NEGF calculations ( Golizadeh-Mojarad and Datta . 200761 ) to SHE experiments all the way 
to room temperature. Another simple way to smear sharp featur es in the loca l spin d en- 



sity due to quantum coherence is averaging over an energy interval (jNikolic et all 12005 d ) , as 



shown in Fig. [91 For more details about different NEGF-based simple models of dephasing in 
quantum transport and their effect on momentum and spin relaxation see the chapter 3 by 
R. Golizadeh-Mojarad and S. Datta in this Volume. 



14 Since current through Biittiker probes is fixed to be zero, then for every electron that enters the lead 
and is absorbed by the corresponding reservoir another one will come with phase memory washed out due to 
equilibration by inelastic effects assumed to take place in the reservoir. Thus, this method introduces inelastic 
effects is a simple phenomenological way. 
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8.5 



Example: Local spin density and spin currents in mesoscopic 
SHE nanostructures 



Figure[9^a) shows NEGF computation (INikolic et all 12005 d ) of spin density within the Rashba 
coupled 2DEG where a ll essential features of the experimentally observed SHE in disordered 
semiconductor devices (IKato et al.l . 12004a ) are obtained but for perfectly clean sample, such 
as: the out-of-plane (S^) component of the spin accumulation develops two peaks of opposite 
signs at the lateral edges of the 2DEG; upon reversing the bias voltage, the edge peaks 
flip their sign; around the left-lead-2DEG interface {S^} is suppressed since unpolarized 
electrons are injected at this contact. By attaching two additional transverse leads at the 
lateral edges in Fig. Mjo), we show how spin-f and spin-j densities will flow through those 
leads in opposite directions to generate spin-Hall current, thereby providing an all-electrical 
semiconductor-based spin injector using ballistic nanostructures with intrinsic SO couplings. 
No discussion of controversial spin currents within the sample is necessary to reach this 
conclusion. Note that simple picture of opposite spin densities flowing in opposite directions 
in Fig. Mjo) is obtaine d in samples smalle r than Lgo, while in larger samples these patterns 
are more complicated (INikolic et all 120061 ) due to multiple spin precession within the sample 
before spin has a chance to exit through the leads. 

The spatial profile of local spin currents corresponding to Fig.[9]^b) is shown in Fig.[9]^c),(d) 
where we separate the integration in Eq. (173bl) for S z bond spin current into two parts 



Ji 



J, 



J, 



S z (neq) 



J mm' 



E F -eV/2 

t f dE 



2 J 2tt 

E F +eV/2 



Trs MG<, m (£)-G< m ,(£))] 



+y / ^Tr s [a,(G<, m (E)-G< m ,(E))] 



E F -eV/2 



plotted as panels (c) and (d), respectively. The states from the band bottom E b to E F — eV/2 
are fully occupied, while states in the energy interval from the electrochemical potential 
Ep — eV/2 (eV > 0) of the right reservoir to the electrochemical potential Ep + eV/2 of the 
left reservoir are partially occupied because of the competition between the left reservoir which 
tries to fill them and the right reservoir which tries to deplete them. The spatial distribution 
of the microscopic spin currents in Fig. M^c) is akin to the vortex-like pattern of bond spi n 
currents within the device that would exist in equilibrium eV p = const. (INikolic et al.l . I2006T ) . 
These do not transport any spin between two points in real space since their sum over any 
cross section is zero. Thus, Fig. [9^d) demonstrates that non-zero spin-Hall flux through the 
transverse cross sections is due to only the wave functions (or Green functions) around the 

Fermi energy. The same sum of (^J^m^^ m (°0 over arbitrary transverse cross section 
(orthogonal to the y-axis) defines the total spin current which, although not the same on 
different cross sections within the sample, flow s into the leads where i t becomes conserved 
and it is measurable [e.g., via the inverse SHE (IHankiewicz et all 12004 )] . 



40 




10 20 30 40 50 60 70 



t so = 0.1t o 






eV=0 






-eV/2 






4 


eV/2 




1 

eV=0 


^8o, 



(c) 



eV/2 



eV=0 



'111' ,t\\, 



r ' 



eV=0 



(d) 



-eV/2 



eV/2 



8a 



eV=0 



' 

>...s!xs..i.:... 

> V \ 1 ■ ' • 

> \ \ \ \ ' 
1 1 \ 1 1 1 

' I I I I 
' I I 1 I 

''111 

' I I I I 



eV=0 



eV/2 



8a 



Figure 9: (a) The out-of-plane component (S^) of the nonequilibrium spin accumulation 
induced by nonlinear (eV = 0.4£q) quantum transport of unpolarized charge current injected 
from the left lead into a two-terminal clean 2DEG nanostructure (of size L = 30a > L^o, a ~ 3 
nm) with the Rashba SO coupling tso = O.lto and spin precession length Lso ~ 15.7a. (b) 
Lateral spin-f and spin- J, densities will propagate in opposite directions through the attached 
transverse ideal (t so = 0) electrodes to yield a linear response spin-Hall current Iy z flowing 
out of 2DEG (L = 8a < L so ), which changes sign upon reversing the bias voltage, (c), (d) 
The spatial distribution of local spin currents in ballistic four-terminal 2DEG bridges with 
the Rashba SO coupling tso = O.lto setting Lso ~ 15.7a and linear response bias voltage 
eV = 10~ 3 to- The local spin current is the sum of equilibrium (persistent) spin current in 
(c), carried by the fully occupied states from — 4t to E F — eV/2, and the nonequilibrium 
(transport) spin current in (d) carried by the partially occupied states around the Fermi energy 
from [ir = Ep — eV/2 (electrochemical potential of the right reservoir) to hl — Ep + eV/2 
(electrochemical potential of the left reservoir). 
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9 COMPUTATIONAL ALGORITHMS FOR REAL(g)SPIN 
SPACE NEGFS IN MULTITERMINAL DEVICES 



While all formulas in Sec. 18.41 for the core NEGF quantities, G(E) and G < (E), can be 
implemented by brute force operations on full matrices, this is typically restricted to small 
lattices of few thousands sites due to computational complexity of matrix operations. For 
example, for a system of size L in d dimensions, the computing time scales as L 3d while the 
memory needed scales as L 2d . By separating system into slices described by much smaller 
Hamiltonian mat rices, and by using the recursive Green function algorithm in serial or parallel 
implementation (jDrouvelis et. all 120061 ) . the complexity can be reduced drastically to L 3d ~ 2 
scaling of the required co mputing time. 



It is often considered (IKazymyrenko and Waintall . 120081 ) that the recursive algorithm can 
be applied only to two-terminal devices and for the computation of total transport quantities 
in the leads. Some alternative algorithms with reduced computational complexity tailore d 
for multiterminal structures have been proposed recently (IKazymyrenko and Waintall . 120081 ) . 
and the rec ursive algori t hm has also been extended to get local charge qu antities within 
the sample (ICresti et all 120031 ; iMetalidis and Brunol . 120051 ; iLassl et all 120071 ) . We discuss in 
Sec. 19.21 how to extend the recursive scheme to four-terminal mesoscopic spin-Hall bridges in 
Fig. [7] that makes possible computation of the spin-Hall conductance on large lattice sizes 
shown in Fig. |8](a). In addit ion, Sec. [9T31 s hows how to obtain local spin densities within 
large two-terminal structures ( Nikolic et al. , 2005 d ) using the recursive-type approach. Since 
the starting point of any algorithm to obtain G(E) and G K (E) for open finite-size system is 
computation of self-energies for the infinite part (ensuring continuous spectrum and dissipa- 
tion) of the system in the form of semi-infinite electrodes attached to the finite-size sample, 
we briefly review recent developments in numerical techniques used to compute the surface 
Green function of semi-infinite leads that generates matrices for "S(E) and 12 < (E). 



9.1 Numerical algorithms for computing self-energy matrices 

Several different numerical techniques are availabl e to evaluate single-part icle Green function 
of a semi-infinite electro de using a lo calized basis (jVelev and Butlerl . 120041 ). For example, the 
so-called Ando method (lAndol . Il99ll ) computes the surface Green function of a homogeneous 



lead, consisting of repeating supercell described by the Hamiltonian H^j where supercells 
are coupled by the Hamiltonian H^j-i, from its exact Bloch propagating modes at energy 
E. This m ethod offers high precis i on numerica l evalu ation of the s elf-energy and can be 
generalized ( Khomyakov et al. . 2005 : Rocha et al. . 2006 ) or accelerated ( Sorensen et all 2008 ) 
to complicated homogeneous periodic systems where Hj j_! is not invertible (required in the 
original Ando algorithm) due to atomistic structure of the leads or more complicated than 
squar e tight-binding lattice [such as the honeycomb one (IChen et all 120071 : IZarbo and Nikolid . 
20071 )]. 

An alternative and approximate method, applicable also to inhomogeneous systems (when 
the lead does not consist of repeating supercells with identic al H, ; , ; ), is the re c ursive (or 
continued fraction) algorithm for the surface Green function (jVelev and Butlerl . 120041 ). It 
connects Green function of a given layer to Green functions of neighboring layers, so that 
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Figure 10: The conductance quantization (a) in the Rashba spin-split infinite clean quantum 
wire modeled on the lattice with three sites per cross section. The wire is decomposed into 
the finite-size central sample attached to two semi-infinite leads where the Rashba hopping 
^so — exists in all three segments. The dashed line in (a) corresponds to quantized 
conductance of the same wire in the absence of the SO coupling tso — 0. Panel (b) plots the 
Rashba-split subband dispersion responsible for this unusual conductance quantization. 



starting from the surface layer and repeating this process until the effect of all other layers on 
the surface layer is taken account yields the surface Green function (in practice, the recursive 
method is executed to finite depth into the lead where the number of considered neighboring 
layers affecting the surface layer is set by the convergence criterion). 

The A ndo m e thod can be executed analytically for simple square lattice lead with no SO 
coupling (IDattal . Il995l ). For comparing mesoscopic SHE calculations to the bulk SHE de- 
scription, it is also useful to have a system where SO-coupled leads are attached to identical 
SO-coupled sample to form an infinite spin-split wire, while two additional ideal (w ith no spin 
interactions) leads are attached in the transverse direction (jSheng and Tingi . 120061 ) . The self- 
energies for semi-infinite electrodes with the Rashba SO coupling can be obtained via the Ando 
method, with doubled size matrices and H^i-i to include spin. We employ this method in 
Fig. EH to provide a simple example of the conductance quantization through three-channel 
infinite spin-split quantum wire that can be used as a building block of the four-terminal 
calculations. The Rashba-split subbands define new conducting cha nnels whose spin polar- 
izatio n properties can be highly nontrivial for strong SO coupling (IGovernale and Ziilickd . 
20041 ). as discussed in Sec. 16.2.11 This is in contrast to standard wires where conducting 



chan nels are always sepa r able st ates |n) <8> |cr), which simplifies the analysis of spin injec- 
tion (jNikolic and Soumal . 12005 a ) while also introduci ng the scattering at the lead-sample 
interf ace which can reduce the spin-Hall conductance (jNikolic et all 12005 a ISheng and Tingi . 



2006). 
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Figure 11: Illustration of the recursive Green function computational algorithm for: a) two- 
terminal setup; and b) four-terminal setup defined on the square tight-binding lattice. 



Recursive Green function algorithm for total spin and charge 
currents in multiterminal nanostructures 



9.2 



For many transport calculations the full retarded Green function Eq. ( 1821) is not required. 
For instance, to compute the transmission matrix ( )4*9l) from lead q to lead p, one only needs 
the Green function matrix elements (G) mm / iCrcr / where the sites m belong to the lead p edge 
supercell and the sit es m' belong to the fi rst supercell of the lead q. We first discuss the 
recursive algorithm (jDrouvelis et. all 120061 ) that evaluates only this portion Gat + i i0 of the 
full retarded Green function matrix for the two-probe setup in Fig. [TTl a). The sample is 
divided into iV layers, and to account for the effect of the leads, their self-energies are added 
onto the layers and N + l. These two layers are the first principal layers of the two leads. 
It is important to emphasize that adding the self-energies directly onto the sample in more 
complicated cases than simple square lattice can lead to incorrect result for the transmission 
eigenvalues. 

The recursive scheme starts by adding the self-energy of the left lead X l onto its principal 
layer directly connected to the sample. The Hamiltonian of any of the N + 2 blocks is labeled 
by h n . The retarded Green function matrix of the n-th isolated block is g n = (E — h n — irj)^ 1 . 
The hopping matrix connecting blocks n and n + 1 is 



-rrhop 

n n,n+l 





H n +i,n 



H n ,n+1 
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We also introduce H n as the matrix which is the sum of the Hamiltonian of n + 1 blocks 
and the lead self-energy added to the 0-th block. The corresponding retarded Green function 
matrix is G n = (E — H n + irf)~ l . One can add the next isolated block h n+1 to H n and obtain 

Assuming that the Green function of n blocks is known, we add the n + 1-th block to 
H n without turning on the coupling H^ +1 between them. The new Hamiltonian, is block 
diagonal 

' Hn 
h n+ i 



H 



(o) 

n+l 



9) 



whose inversion leads to 



G 



(o) 

n+l 



G„ 






ln+1 



(90) 



To find the Green function G n+ i it is necessary to add the matrix elements connecting the 
n + 1-th block to H n 

H n +i )?1 h n+ i 

The Green function of n + 2 blocks can be found by using the Dyson equation 



H 



n+l 



(91) 



Gf~i{0) ,-,■■>> XT Hop r , 

n+l — W+l "+" ^n+l n n,n+l^n+l' 



(0) Trhop 



(92) 



All matrices in this equation are of the dimension n + 2 in the block coordinate space. Equa- 
tion ( |92i) can be used to obtain the recurrence relationships for the Green function blocks 
we are interested in. To get the conductance of the two-terminal device we need to compute 
(N + l | G jv+i | 0) = Gat + i(A^ + 1, 0) which connects the right and the left lead. 

In the multiterminal case, we choose to start from the Hamiltonian of the block of sites in 
the center of the sample and then add one-by-one additional blocks of the same shape until 
we reaches the attached electrodes, as sketched in Fig. [TlTb). Therefore, in this case we need 
the Green function elements GAr + i(iV + 1, N + 1). From Eq. (1921) we find 



G%(v + i.ii I } • = •(// 1 



f~((0) Trhop (~i 
^n+l-^Vn+l^rc+l 



n + l 



(93) 



G n+ i(n + + 

By inserting the identity matrix in the subspace of block coordinates into Eq. ( |93i) 

I = \m) (m\ , (94) 



the second term in Eq. ( 1931) becomes 



71 + 1 

It is easy to see from Eqs. 



GitiEN H < n p +1 ^ |/> (Z|G n+1 



71+1 



and (JHD that 



(n + l|G 



(o) 

n+l 

rhop 



\m) 



8. 



n+l,m&n+l > 



(95) 



(96a) 
(96b) 
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Equation (1961) can be substituted into Eq. (195]) to yield the Dyson Eq. (1921) for G„ + i(n+l, n+1) 

G n+ i(n + l,n + 1) = gn+1 + gn+lH n+ i )n G n+1 (n,n + l), (97) 

in which G n+ i(n, n + 1) is an unknown quantity. By computing the matrix elements of the 
terms in the Dyson Eq. ( 1921) between \n) and \n + 1) we find that 

G n+1 (n,n + 1) = G n \Jli Tl)H n rl +i G n+1 (n + l,n + l). (98) 

Finally, combining Eqs. flDTj) and (19"H1) leads the following recurrence relationship 

G n+1 (n + 1, n + 1) = — % — • (99) 

Here the matrix is the identity matrix in the subspace of the n + 1-th block and the 
denominator is the matrix to be inverted. 

If the recursive scheme in Fig. fTTT b) is used, we need to obtain Gjv +1 (iV + 1, N + 1) 
from which the transmission matrix can be computed. For the two-terminal case, the recur- 
sive scheme sketched in Fig. UTT a) is more efficient computationally. Computing the matrix 
elements of the Dyson Eq. ( l92i) between the states \n + 1) and \m) gives 

G n+ i(n + l,m) = G n+ i(n + 1, n + l)H n+ i in G n (n, m) (100) 

where m < n + 1. 

We summarize the general recursive Green function algorithm for the computation of 
(spin-resolved) transmissions, and therefore spin or charge conductances, in the form of the 
following steps: 

1. Calculate the self-energies for the leads, as described in Section [9~T1 

2. Identify the N + 2 blocks for the recursive method as shown in Fig. [TT] (the number of 
sites of each block might be different). 

3. Initialize the recurrence by computing Go(0, 0) = go, where go = (E — H — iff)' 1 
is the Green function of the first block. This block has the self-energy added to it in 
the two-terminal case. In the multiterminal case, the first block can be chosen to be a 
rectangle of the size iV"i 0) x ^ 0) sites, as illustrated in Fig.EHb). A difficulty may arise 
in the multiterminal cases where the matrix E — ho can turn to be singular for certain 
choices of the shape of the sample or its parameters. 

4. Using the formulas ( 1991) and ( jlQOj) compute the Green functions Gjv + i(iV + 1, JV + 1) 
and G N+1 (N + 1,0). 

5. Obtain the retarded Green function matrix elements G pq connecting each pair of leads 
p and q, either from Gjv + i(iV + 1, 0) in the two-terminal case, or 

G7v + i(iV + 1,N + 1) in the multiterminal case. 

6. Compute the transmission matrix using Eq. 
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9.3 Recursive Green function algorithm for local spin and charge 
densities 

The computation of nonequilibrium spin and charge densities using Eq. flHUj) and (ITUl) . respec- 
tively, requires to obtain G < . This has the highest computational complexity since one has 
to perform the matrix multiplication in Eq. flHTj) . Moreover, this matrix multiplication has 
to be done for all the energy points in Eq. ( IHOi) . Here we introduce a method for recursive 
calculation of the spin and charge nonequilibrium densities within the NEGF formalism which 
makes it possible to greatly reduce this computat ional complexity and, thereby, analyze spin 



densities in sizable samples (jNikolic et all l2005d ) . 



In order to illustrate the essential steps of this method, it is enough to consider a simple 
two-probe setup in Fig. HIT a). The two-probe system considered here is made of iV + 2 blocks, 
where the 0-th block is the first left lead supercell onto which the left lead self-energy is added, 
and N + 1-th block is the first right lead supercell onto which the right lead self-energy is 
added. The voltage applied onto the left and right lead is Vl and V R , respectively. The lesser 
self-energy matrix in Eq. flHTT) can be rewritten as 



£<(£) = iT L {E-eV L )f{E-eV L )+iT R {E-eV R )f{E-eV R ) 
iT L (E-eV L )f(E-eV L ) ... 

... iT R (E-eV R )f(E-eV R ) 



(101) 



where the only non-zero elements are the first and the last on the diagonal. This simplifies 
greatly the matrix multiplication in Eq. (1ST]) , since only multiplication by the non-zero matrix 
elements of S < is required 

£<(0,0) = iT L (E-eV L )f(E-eV L ), (102a) 
S<(JV + l,iV + l) = iT R (E-eV R )f(E-eV R ). (102b) 

In order to calculate the nonequilibrium spin density Eq. (IHUI) or charge density Eq. (ITU]) , only 
the diagonal elements of G < (E) are required. Using Eqs. (1101)) and (11021) in Eq. (IHTj) . leads 
to 

G<(m,m) = G(m,N + 1)E < (JV+ l,N+ l)G f (iV + l,m) 

+ G(m,0)S < (0,0)G t (0,m), (103) 

where m labels the coordinates of the m-th block. 

The retarded Green function matrices G(m, N + 1) and G(m, 0) can be found using the 
Dyson Eq. (I9"2"j) introduced in Sec. 19.21 To find G(m, N + 1), one starts the recurrence from 
the left lead and finds 

G n+ i(m, n + 1) = G n (m, n)H n>n+ iG n+ i(n + 1, n + 1), (104) 

where m < n. The Green function block G n+ i(n + l,n+ 1) can be obtained from Eq. (l9TJj) . 
To find G(m, 0), the recurrence has to start from the right lead and we get 

G n - 1 (m,n - 1) = G n (m,n)H niri _iG n _i(72 - l,n - 1), (105) 
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where m > n. The system of equations is finally closed by using 

Gn - l(n -^ n -^ = g -i h In C(nnM ' (1 ° 6) 

In this method, one has to compute recursively the retarded Green functions G(m, N + l) and 
G(m, 0) for all m = 0, N + 1. It is not as computationally costly as the full Green function 
matrix inversion, and it can be extended to compute bond spin and charge currents. In that 
case, we also have to find G(m ± 1, N + 1) and G(m ±1,0). 



10 CONCLUDING REMARKS 

The pure spin currents, which arise when equal number of spin-f and spin- j electrons propa- 
gate in opposite directions so that net charge current carried by them is zero, have recently 
emerged as one of the major topics of both metal and semiconductor-based spintronics. This 
is due to the fact that they offer new realm to explore fundamental spin transport phenomena 
in solids, as well as to construct new generation of spintronic devices with greatly reduced 
power dissipation and new functionality in transporting either classical or quantum (such as 
flying spin qubits) information. Their generation (which is currently experimentally achieved 
mostly through spin pumping, non-local spin injection using lateral spin valves, and SHE), 
control, and demanding detection typically involves quantum-coherent spin dynamics. In 
particular, the recently discovered SHE and the inverse SHE have offered some of the most 
efficient scheme to generate and detect pure spin currents, respectively. However, they have 
also posed numerous challenges for their theoretical description and device modeling since 
conventionally defined spin current operator does not satisfy standard continuity equation for 
spin when intrinsic SO couplings (which act as homogeneous internal momentum-dependent 
magnetic fields affecting the band structure and causing spin precession) are present. Even 
with modified spin current definitions, this fundamental problem makes difficult easy con- 
nection between theoretically computed spin conductivities of infinite homogeneous systems 
and usually experimentally measured spin densities for devices with arbitrary boundaries and 
attached external electrodes. An attempt to give up on the theoretical description of spin 
flow through spin current altogether, by using spin density within the diffusion equation for- 
malism (which is expected to provide a complete and physically meaningful description of 
spin and charge diffusive transport in terms of position and time dependent spin and charge 
densities), leads to another set of obstacles. That is, the explicit solution of the diffusion 
equation strongly depends on the boundary conditions (extracted for charge diffusion from 
conservation laws) where spin non- conservation, manifested as ballistic spin precession at the 
edges, cannot be unambiguously matched with the diffusive dynamics in the bulk captured 
by the diffusion equation for length scales much longer than the mean free path. 

This chapter of the Handbook provides a summary of such fundamental problems in spin 
transport through SO-coupled semiconductors and possible resolution through spin-dependent 
NEGF formalism that can handle both ballistic and diffusive regimes, as well as phase-coherent 
and dephasing effects, in realistic multiterminal nanostructures operating at different mobil- 
ities and temperatures. The central quantities of this approach — total spin currents flowing 
out of the device through ideal (interaction free) electrodes, local spin densities within the 
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sample and along the edges, and local spin currents (whose sums within the leads gives total 
spin currents) — can be used to understand expe riments or model novel spintronics devices. 



2004) for comput a tiona l 



20081 : iRoth et all |2009j) 



For example, NEGF approach is indispensable (IHankiewicz et al. 
modeling of generator-detector experimental setups (jBriine et al 
where direct and inverse SHE multiterminal bridges are joined into a single circuit to sense 
the pure spin-Hall current through voltages induced by its flow through regions with SO inter- 
actions. The chapter also provides extensive coverage of relevant technical and computational 
details, such as: the construction of retarded and lesser nonequilibrium Green functions for 
SO-coupled nanostructures attached to many electrodes; computation of relevant self-energies 
introduced by the lead-sample interaction; and accelerated algorithms which make possible 
spin transport modeling in device whose size is comparable to the spin precession length of 
few hundreds of nm where SHE response to injected unpolarize d charge current is expected 



to be optimal (INikolic et all 12005 a ; iMoca and Marinescul . 120071 ) . 
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